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CONVERSION  FACTORS,  U.S.  CUSTOMARY  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 


U.S.  customary  units  of  measurement  used  in  this  report  can  bo  converted 
to  notric  (SI)  units  as  follows: 


Multiply 

by 

To  obtain 

inches 

25.4 

millimeters 

2.S4 

centimeters 

square  inches 

6.452 

square  centimeters 

cubic  inches 

16.39 

cubic  centimeters 

feet 

30.48 

centimeters 

0.3048 

meters 

square  feet 

0.0929 

square  meters 

cubic  feet 

0.0283 

cubic  meters 

yards 

0.9144 

met  ars 

square  yards 

0.836 

square  meters 

cubic  yards 

0.7646 

cubic  meters 

miles 

1.6093 

kilometers 

square  miles 

259.0 

hectares 

Vncts 

1.8532 

kilometers  per  hour 

acres 

0.4047 

hectares 

foct-pounds 

1.3558 

newton  meters 

millibars 

1.0197  x 

10"  3 

kilograms  per  square  centimeter 

ounces 

28.35 

grams 

pounds 

453.6 

grams 

0.4536 

kilograms 

ton,  long 

1.0160 

metric  tons 

ton,  short 

0.9072 

metric  tons 

degrees  (angle) 

0.1745 

radians 

Fahrenheit  degrees 

5/9 

Celsius  degrees  or  Kelvins1 

To  obtain  Celsius  (C)  temperatui'c  readings  fro r Fahrenheit  (F)  readings, 
use  formula:  C = (5/9)  (F  -32). 

To  obtain  Kelvin  (K)  readings,  use  formula:  K = (5/9)  (F  -32)  + 273.15. 
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SYMBOLS  AND  DEFINITIONS 


A 

ANG 

CD 


CL 

clear 

s, 

d 

Dia 

F 

FD 

(vh 

(Vv 

Fh 

w 

FI 

(FI5h 

(FIJV 


Fi 

Foh<V 

Fov«V 


F (0.) 
v 1 


projected  area  of  pipe  section 
orientation  angle  with  respect,  to  wave  crests 
coefficient  of  drag 
coefficient  of  lift 

coefficient  of  transverse  force  due  to  eddy  shedding 

bottom  clearance 

coefficient  of  mass 

Stillwater  depth 

pipe  diameter 

total  wave-induced  force 

drag  force 

horizontal  component  of  drag  force 

vertical  component  of  drag  force 

iiorizontal  component  of  total  wave  force 

calculated  horizontal  force  at  position  0.  in  wave  cycle 

inertial  force 

horizontal  component  of  inertial  force 
vertical  component  of  inertial  force 
lift  force 

transverse  "lift"  force  due  to  eddy  shedding 
observed  horizontal  force  at  position  0.^  in  wave  cycle 
observed  vertical  force  at  position  0^  in  wave  cycle 
vertical  component  of  total  wave  force 
calculated  vertical  force  at  position  0^  in  wave  cycle 
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SYMBOLS  AND  DEFIN  TTQNS~ Continued 


H 

k 

L 

T 

t 

u 

u 

max 

V 
v 

v 

max 

z 

3u/'ot 

3v/3t 

0 

V 

P 

4> 

Input 

ANGLE 


wave  height 

negative  fraction  of  lift  force  cycle 

wavelength 

wave  period 

time  since  last  wave  crest  passed  over  center  of  pipe 
section 

horizontal  component  of  water  particle  velocity  if  pipeline 
was  absent 

maximum  horizontal  water  particle  velocity  if  pipeline  was 
absent 

volume  of  fluid  displaced  by  pipe  section 

vertical  component  rc  water  particle  velocity  if  pipeline 
was  absent 

maximum  vertical  water  particle  velocity  if  pipeline  was 
absent 

vertical  distance  of  center  of  pipe  section  above  oottom 

horizontal  component  of  water  particle  acceleration  if  pipe- 
line was  absent 

vertical  component  of  water  particle  acceleration  if  pipe- 
line was  absent 

2irt/T  = position  of  wave  cycle  over  center  of  pipe  section 
with  respect  to  time 

kinematic  viscosity  of  fluid 

mass  density  of  fluid 

phase  shift  of  maximum  lift  forces  with  respect  to  wave  cycle 
Computer  Programs 


Parameters : 

orientation  angle 

calibration  factor  for  manual  digitizer 
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SYMBOLS  AND  DEFINITIONS— Continued 
downward  force  calibration  factor 
upward  force  calibration  factor 
bottom  clearance 

downward  force  calibration  factor 
pipe  diameter 

negative  wave  (trough)  calibration  factor 

wave  force  readings 

zero  point  of  wave  force  record 

wave  surface  readings 

number  of  wave  force  readings 

wave  period 

upward  force  calibration  factor 
positive  wave  (crest)  calibration  factor 
zero  point  of  wave  record 
length  of  pipe  test  section 


XF  amplification  factor  for  force  record 

XW  amplification  factor  for  wave  record 

YI(I)  wave  surface  readings 

Program  Variables: 

ANG  orientation  angle  (in  radians) 

ANGLE  orientation  angle  (in  degrees) 

CDH  horizontal  coefficient  of  drag 

CDV  vertical  coefficient  of  drag 

CL  bottom  clearance 


jit  II.?!,', 
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SYMBOLS  AND  DEFINI TIONS — Continued 


II  \ 

I 


coefficient  of  lift  (calculated  using  horizontal  velocity 
in  direction  of  wave  advance  end  projected  area  in  plane 
parallel  to  the  pipeline  axis) 

coefficient  of  lift  (calculated  using  horizontal  velocity 
in  direction  of  wave  advance  and  projected  area  in  plane 
normal  to  the  direction  of  wave  advance) 

coefficient  of  lift  (calculated  using  the  component  of  the 
horizontal  velocity  in  the  direction  perpendicular  to  the 
pipeline  axis  and  the  projected  area  in  the  plane  parallel 
to  the  pipeline  axis) 

horizontal  coefficient  of  mass 

vertical  coefficient  of  mass 


Stillwater  depth 

pipe  diameter 

1/2  p A u2 

max 

1/2  P A v2 

' me  v 


FH(1)  calculated  horizontal  wave  force 

FI (I)  measured  wave  force  readings  (in  grams  for  two-dimensional 

data;  in  10-grams  for  three-dimensional  data) 


FLV  1/2  p A u2 

K max 

FMAX  maximum  positive  wave  force  (measured) 

FMH  P V (3u/3t)max 

FMIN  maximum  negative  wave  force  (measured) 

FMV  p V (3v/0t) 

FP(1)  measured  wave  force  readings  (in  pouuds) 

«r 

FV(I)  calculated  vertical  wave  force 


HI  (I) 


wave  height 

wave  surface  profile  readings 
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SYMBOLS  AND  DEFINIi  IONS— Continued 


PHI  phase-shift  parameter  <J>  of  modified  lift  force  equation 

PI  TT 

R mass  density  of  water 

RES (I)  difference  between  measured  wave  force  and  calculated  wave 
force 


SF  wave  force  averaged  through  wave  cycle 

T wave  period 

U maximum  horizontal  water  particle  velocity 

XC  length  of  pipe  section 

XK  parameter  K of  modified  lift  force  equation 

XL  wavelength 

Z V vertical  distance  from  bottom  to  center  of  pipe  section 

Tabulated  Experimental  Data 

ANG  orientation  angle  of  pipeline  with  respect  to  wave  crests 

CDH  horizontal  coefficient  of  drag 

CDV  vertical  coefficient  of  drag 

CLER  bottom  clearance 

CLV  coefficient  of  lift  (calculated  using  horizontal  velocity 

in  direction  of  wave  advance  and  projected  area  in  plane 
parallel  to  the  pipeline  axis) 

CLVA  coefficient  of  lift  (calculated  using  horizontal  velocity 

in  direction  of  wave  advance  and  projected  area  i;i  the 
plane  normal  to  the  direction  of  wave  advance) 

CLVU  coefficient  of  lift  (calculated  using  the  component  of  the 

horizontal  velocity  in  the  direction  perpendicular  to  the 
pipeline  axis  and  the  projected  area  in  the  plane  parallel 
to  the  pipeline  axis) 

CMH  horizontal  coefficient  of  mass 
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CMV  vertical  coefficient  of  mass 

DIA  pipe  diameter 

FAVC  average  horizontal  force  (averaged  over  complete  wave  cycle) 

H wave  height 

K.  parameter  k of  modified  lift  force  equation 

L wavelength 

PHT  phase  shift  parameter  <j)  of  modified  lift  force  equation 

T wave  period 

UMAX  Maximum  horizontal  component  of  water  particle  velocity  at 

center  of  pipe  section  if  absent 
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FORCES  EXERTED  BY  WAVES  ON  A PIPELINE 
AT  OR  NEAR  THE  OCEAN  BOTTOM 

by 

George  L.  Bowie 
I.  WAVE  FORCE  ANALYSIS 


1 . Wave  Force  Components  on  Pipelines  Near  the  Bottom. 

The  most  common  method  of  analyzing  wave  forces  on  pipelines  is 
the  application  of  the  Morison  equation  (Morison,  et  al.,  1950).  Using 
this  approach,  the  total  wave- induced  force  on  a pipeline  can  be  broken 
into  several  components,  depending  on  whether  the  components  are  due  to 
the  water  particle  velocities  or  accelerations.  These  force  components 
can,  in  turn,  be  separated  into  horizontal  and  vertical  components 
by  using  the  horizontal  and  vertical  components  of  the  water  particle 
velocities  and  accelerations  in  their  respective  force  equations. 

Where  there  is  no  lift  effect  and  no  eddy-induced  forces,  the  vertical 
component,  Fv,  of  the  total  wave  force  is 


Fv  ■ (Fi)y  - (FD)v  = <4,  p V + 1/2  CD  p A v|v|  (1) 

and  the  horizontal  component,  Fh,  is 

Fh  * (Fl)h  + (FD)h  - CM  p V |H-  + 1/2  CD  P A u|u|,  (2) 


where 

Wv 

(Fl)h 

(pD)v 

Wh 


u 


3v 

3t 


vertical  component  of  inertial  force 

horizontal  component  of  inertial  force 

vertical  component  of  drag  force 

horizontal  component  of  drag  fcrce 

vertical  component  of  water  particle  velocity  if 
pipeline  was  absent 

horizontal  component  of  the  water  particle  velocity 
if  pipeline  was  absent 

vertical  component  of  water  particle  acceleration 
if  pipeline  was  absent 


3u 

3t 

A 

V 

P 

Si 


horizontal  component  of  water  particle  acceleration 
if  pipeline  was  absent 

projected  area  of  pipe  section 

volume  of  fluid  displaced  by  pipe  section 

mass  density  of  fluid 

coefficient  of  mass 

coefficient  of  drag 


For  a pipeline  located  near  the  ocean  bottom,  the  water  particle 
orbits  are  flattened  parallel  to  the  boundary.  Assuming  a horizontal 
bottom,  the  vertical  motions  of  the  water  particles  are  small  in  com- 
parison to  the  horizontal  motions,  especially  in  shallow-water  depths 
relative  to  the  wavelength.  As  a result,  the  vertical  components  of 
the  water  particle  velocities  and  accelerations  are  much  smaller  than 
the  horizontal  components,  and  correspondingly  the  vertical  components 
of  the  drag  and  inertial  forces  will  be  smaller  than  the  analogous 
horizontal  forces. 


Since  the  water  particles  at  the  bottom  are  effectively  oscillating 
in  a horizontal  plane,  the  vertical  excursions  of  the  water  particles 
will  generally  be  less  than  the  diameter  of  a submarine  pipeline  lying 
on  or  near  the  bottom.  Therefore,  the  vertical  drag  forces  are 
generally  insignificant,  and  could  probably  be  neglected  from  the 
vertical  wave  force  equation. 

Pipelines  near  the  bottom  are  subject  to  vertical  lift  forces. 

These  forces  are  the  result  of  the  asymmetric  distortion  of  the  flow 
field  due  to  the  proximity  of  the  bottom  boundary,  which  induces  dif- 
ferences in  the  horizontal  flow  velocities  and  corresponding  pressure 
distribution  over  the  top  and  bottom  of  the  pipeline.  Since  the  water 
particle  velocities  near  the  bottom  are  at  a maximum  in  the  horizontal 
plane,  the  lift  forces  induced  by  these  horizontal  motions  will  gener- 
ally be  the  predominant  force  acting  in  the  vertical  direction. 

Transverse  "lift"  forces  due  to  eddy  shedding  may  also  be  an 
important  component  of  the  vertical  wave  force,  since  these  forces  are 
also  due  to  the  horizontal  water  particle  velocities  and  excursions 
which  are  maximum  in  the  horizontal  direction.  Certain  values  of  the 
Keulegan-Carpenter  parameter  and  Reynolds  number  must  be  attained  for 
the  eddy  release  phenomenon  to  occur.  The  proximity  of  the  bottom 
boundary  will  probably  have  some  effect  on  the  formation  and  release 
of  the  eddies,  both  because  it  is  a solid  boundary,  and  because  it 
affects  the  orbital  motions  of  the  water  particles  induced  by  the 
wave  action. 
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Although  the  eddy-induced  component  of  the  vertical  wave  force  may 
be  significant  wnen  compared  to  the  relatively  small  vertical  drag  and 
inertial  forces,  the  experimental  results  of  this  investigation  show 
that  the  eddy-induced  lift  forces  are  much  smaller  than  the  "Bernoulli- 
type”  lift  forces  for  pipelines  located  near  the  bottom.  At  large 
clearances  above  the  bottom  where  the  Bernoulli-type  lift  effect  becomes 
negligible,  the  transverse  lift  forces  due  to  eddy  shedding  may  become 
a significant  component  of  the  total  vertical  force.  At  the  same  time, 
as  the  pipeline  is  raised  farther  from  the  bottom  boundary,  the  verti- 
cal inertial  and  drag  forces  also  become  more  significant. 

The  vertical  component  of  the  total  wave-induced  force  acting  on  a 
pipeline  near  the  ocean  bottom  thus  consists  of  four  components — the 
lift  force,  the  inertial  force,  th.  drag  force,  and  the  transverse  lift 
force  due  to  eddy  shedding.  Using  the  Morison  approach,  the  total  ver- 
tical wave  force  is  expressed  as  the  sum  of  these  components : 

Fv  - FL  * <FI>V  * <FD>v  * V <3> 

4 

where  is  the  lift  force  and  F^’  is  the  transverse  lift  force  due  to 
eddy  shedding . 

2.  Wave- Induced  Lift  Forces. 

Cansider  a pipeline  in  contact  with  e.  horizontal  rigid,  impervious 
bottom.  Water  cannot  flow  between  the  pipe  and  the  bottom  boundary,  so 
the  flow  must  be  diverted  over  the  top  of  the  pipe.  The  asymmetrical 
distortion  of  the  flow  field  results  in  maximum  velocities  over  the  top 
of  the  pipe  section  and  minimum  velocities  oyer  the  bottom,  with  zero 
velocities  at  the  stagnation  point  on  the  upstream  side  of  the  pipe 
bottom  at  the  point  of  contact  with  the  sea  floor.  Correspondingly, 
the  associated  pressure  distribution  will  induce  an  upward  lift  force 
for  any  velocity  field  acting  on  the  pipeline.  The  stagnation  pressure 
at  the  bottom  of  the  pipe  section  will  increase  with  increasing  veloc- 
ity, while  simultaneously  the  pressure  distribution  over  the  top  of  the 
pipeline  will  decrease  with  the  increased  velocities  of  the  flow  di- 
verted over  the  top  of  the  pipe  section.  The  wave-induced  lift  forces 
will  thus  act  in  the  upward  direction  throughout  the  wave  cycle,  in- 
creasing with  the  horizontal  water  particle  velocities  to  maximum  mag- 
nitudes under  the  crests  and  troughs  of  the  passing  waves,  and 
diminishing  to  zero  at  the  points  of  horizontal  flow  reversal. 

In  contrast,  a pipeline  located  at  a small  clearance  above  the 
bottom  boundary  is  subject  to  a more  complex  type  of  lift  phenomenon. 

At  the  phase  in  the  wave  cycle  where  the  horizontal  component  of  ^e 
water  particle  velocity  reverses  direction,  the  horizontal  velocity 
over  the  pipeline'  is  approximately  zero.  As  the  wave  crest  or  trough 
begins  to  approach  the  pipeline,  the  wave-induced  horizontal  velocities 
are  initially  low,  inducing  unrestricted  flow  at  low  velocities  over 
both  the  top  and  bottom  of  the  pipeline.  However,  the  water  flows 
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faster  through  the  bottom  clearance  constriction  than  over  the  top  of 
the  pipeline,  so  the  corresponding  differences  in  the  pressure  distri- 
bution exert  a downward  (negative  lift)  force  toward  the  bottom  bound- 
ary (Fig.  1,  a). 

At  first,  the  negative  lift  force  will  increase  with  the  increas- 
ing horizontal  water  particle  velocities  of  uhe  approaching  wave,  since 
the  flow  velocities  increase  at  a faster  rate  through  the  bottom  clear- 
ance constriction  than  over  the  top  of  the  pipeline,  thus  producing 
larger  differences  in  the  corresponding  pressure  distributions  over  the 
top  anu  bottom  of  the  pipe  section  (Fig.  1,  b) . 

This  continues  until  viscous  effects  begin  to  restrict  the  flow 
through  the  narrow  bottom  clearance.  For  a given  small  clearance  and 
a given  amount  of  energy  in  the  horizontal  water  particle  velocities 
approaching  tne  pipeline,  the  velocities  and  flow  rates  of  a viscous 
fluid  through  the  bottom  clearance  constriction  can  attain  only  certain 
maximum  values.  Thus,  a "choking”  effect  is  exerted  on  the  restricted 
flow  through  the  small  bottom  clearance,  and  the  remainder  of  the  wave- 
induced  flow  is  forced  to  flow  over  the  top  of  the  pipe  section.  Cor- 
respondingly, the  stagnation  point  will  shift  downward,  increasing  the 
pressure  on  the  lower  upstream  side  of  the  pipeline.  The  larger  the 
proportion  of  the  flow  diverted  over  the  top  of  the  pipe,  the  lower  the 
stagnation  point. 

At  the  same  time,  the  increasing  velocities  associated  with  the 
approaching  wave  crest  cause  the  restricted  flow  through  the  bottom 
clearance  to  form  a turbulent  jet  with  the  generation  of  eddies  behind 
the  jet.  The  generation  of  increased  turbulence  and  eddies  results  in 
an  energy  loss  in  the  water  flowing  through  the  bottom  constriction, 
decreasing  the  velocities  under  the  pipe  section  behind  the  jet. 

The  above  effects  associated  with  the  choking  phenomenon  limit  the 
maximum  flow  velocities  and  minimum  pressures  under  the  bottom  side  of 
the  pipe  section.  In  contrast,  the  unrestricted  flow  velocities  over 
the  top  of  the  pipeline  increase  freely  with  the  increasing  horizontal 
velocities  of  the  advancing  wave.  The  increased  part  of  the  approach- 
ing flow  that  is  diverted  over  the  top  of  the  pipe  section  due  to  the 
shift  in  stagnation  point  produces  a further  increase  in  the  flow 
velocities  over  the  top.  Correspondingly,  the  pressure  distribution 
over  the  top  side  of  the  pipeline  decreases  at  a faster  rate  than  the 
associated  pressures  along  the  bottom  side,  so  the  negative  lift  force 
gradually  decreases  and  eventually  becomes  positive  (Fig.  1,  c,  d,  and 
e) . 


At  this  stage,  the  upward  lift  force  becomes  larger  as  the  hori- 
zontal velocities  acting  on  the  pipeline  increase  further  with  the 
advancing  wave  crest  or  trough  (Fig.  1,  f). 


As  the  wave  crest  or  trough  passes,  this  series  of  steps  in  the 
lift  force  phenomenon  is  reversed.  The  horizontal  velocities  approach- 
ing the  pipe  section  begin  to  decrease,  resulting  in  a decrease  in  the 
positive'  lift  force  exerted  on  the  pipelin.:.  As  the  velocities  decrease 
further  with  the  passing  wave,  the  flow  under  the  pipe  section  begins 
to  i econe  less  res^^cied.  The  choking  effect  thus  decreasos,  and  the 
turbulence  and  eddies  near  the  bottom  clearance  gradually  diminish.  As 
the  flow  wider  the  pipe  section  ceases  to  be  restricted,  less  of  the 
horizontal  flow  approaching  the  pipeline  is  forced  to  flow  over  the  top 
of  the  pipe,  so  the  stagnation  point  will  accordingly  shift  upward, 
closer  to  the  center  of  the  pipe  section. 

The  flow  velocities  decrease  simultaneously  over  the  top  and  bottom 
of  the  pipeline  as  the  wave  passes,  but  the  rate  of  decrease  is  faster 
over  the  top  of  the  pipe  than  in  the  vicinity  of  the  bottom  constric- 
tion. The  positive  lift  force  decreases  until  eventually,  the  flow 
velocities,  location  of  the  stagnation  point,  and  associated  pressure 
distribution  are  such  that  the  pressure  integrated  over  the  pipe  sec- 
tion again  results  in  a negative  lift  force.  Hie  downward  lift  force 
then  increases  as  the  flow  through  the  bottom  clearance  becomes  less 
restricted  with  the  decreasing  velocities  of  the  passing  wave. 

This  lift  phenomenon,  as  shown  in  Figure  2 for  a passing  wave  crest, 
is  repeated  twice  during  each  wave  cycle  as  the  direction  of  the  wave- 
induced  horizontal  velocities  reverses  under  the  crests  and  troughs  of 
the  passing  waves. 

In  reality,  the  horizontal  flow  reversal  occurs  almost  instanta- 
neously, so  the  negative  lift  force  does  not  return  to  zero  at  the 
point  of  zero  velocity  when  the  flow  reverses  through  the  bottom  clear- 
ance constriction.  The  instant  of  zero  velocity  occurs  only  at  the 
center  of  the  pipe  cross  section  (the  reference  point) . Since  the 
pipeline  has  a finite  diameter,  the  wave- induced  flow  acting  on  the 
pipe  section  at  any  instant  includes  the  sum  of  the  flow  condi- 
tions induced  by  the  part  of  the  wave  covering  the  entire  diameter  of 
the  pipeline.  So  instead  of  going  to  zero  with  the  passing  wave  crest, 
and  then  increasing  initially  with  the  approaching  trough,  the  lift 
force  remains  negative  during  the  period  of  minimal  velocities  as  the 
flow  reverses  under  the  pipe  section. 

In  a similar  manner,  the  lift  force  does  not  become  positive  as 
soon  as  the  choking  effect  occurs  in  the  bottom  clearance  constriction. 
The  development  of  the  choking  phenomenon  involves  the  formation  of  a 
turbulent  jet  through  the  constriction,  and  a downward  shift  in  the 
stagnation  point  as  more  water  is  diverted  over  the  top  of  the  pipe 
with  increasing  restriction  of  the  flow  through  the  clearance.  The 
corresponding  changes  in  the  velocities,  flow  pattern,  and  associated 
pressure  distribution  over  the  top  and  bottom  of  the  pipe  section  pro- 
duce the  transition  from  negative  to' positive  lift.  This  process 
requires  some  small  but  finite  amount  of  time.  Conversely,  the  reversal 
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of  these  processes  with  the  decreasing  velocities  of  the  passing  wave 
crest  also  involves  a small  but  finite  amount  of  time.  Thus,  there  will 
be  a slight  timelag  in  the  point  of  maximum  positive  lift  with  reference 
to  the  instant  of  maximum  velocity  as  the  wave  crest  (or  trough)  passes 
over  the  reference  point.  The  smaller  the  amount  of  positive  lift  rela- 
tive to  the  amount  of  negative  lift,  and  the  later  the  positive  lift 
occurs  in  the  wave  cycle,  the  greater  the  timelag. 

An  example  of  the  lift  force  phenomenon  over  a complete  wave  cycle 
for  a small  bottom  clearance  is  shown  in  Figure  3. 

For  a given  pipe  diameter  and  wave  condition,  as  the  bottom  clearance 
is  increased,  higher  velocities  are  necessary  to  produce  the  choking 
effect  in  which  the  flow  becomes  restricted  through  the  bottom  clearance 
constriction.  Thus,  as  the  bottom  clearance  is  increased,  the  flow  under 
the  pipeline  begins  to  become  restricted  closer  to  the  approaching  wave 
crest  or  trough,  where  the  horizontal  velocities  are  at  a maximum;  this 
choking  effect  also  diminishes  soon  after  the  wave  crest  or  trough  has 
passed.  Therefore,  as  the  bottom  clearance  is  increased,  the  downward 
lift  force  occurs  during  a larger  part  of  the  wave  cycle. 

At  the  same  time,  larger  clearances  permit  greater  maximum  velocities 
and  corresponding  lower  pressures  under  the  pipe  section.  Since  higher 
flow  rates  are  possible  under  the  pipe  section,  less  of  the  wave-induced 
flow  must  be  diverted  over  the  top  of  the  pipeline.  As  a result  of 
these  changes,  the  negative  lift  forces  reach  a greater  magnitude  before 
the  choking  effect  begins,  and  these  maximums  are  attained  later  in  the 
wave  cycle. 

Correspondingly,  the  upward  lift  forces  occur  during  a smaller  part 
of  the  wave  cycle,  and  the  maximum  magnitude  these  forces  attain  decreases 
with  increasing  bottom  clearance.  These  maximum  values  are  also  reached 
later  in  the  wave  cycle. 

If  the  clearance  is  increased  further,  a point  is  eventually  reached 
at  which  the  clearance  is  large  enough  so  that  the  choking  effect  does 
not  occur.  At  this  stage,  the  velocities  are  higher  through  the  bottom 
clearance  constriction  than  over  the  top  of  the  pipeline  during  the 
entire  wave  cycle.  So  the  associated  pressure  distribution  results  in 
a negative  lift  force  throughout  the  wave  cycle,  with  maximum  downward 
forces  occurring  under  the  crests  and  troughs  of  the  passing  waves.  The 
negative  lift  diminishes  to  zero  at  the  points  of  horizontal  flow  reversal. 

As  the  bottom  clearance  is  increased  further,  the  downward  lift  effect 
is  gradually  reduced.  The  phase  of  the  force  cycle  relative  to  the  wave 
cycle  remains  the  same,  but  the  magnitude  decreases.  Eventually,  a 
point  is  reached  where  the  bottom  clearance  no  longer  acts  as  a constric- 
tion to  the  wave-induced  flow.  The  flow  pattern  becomes  approximately 
symmetrical,  and  the  increased  velocities  of  the  horizontal  flow  diverted 
over  the  top  and  bottom  of  the  pipeline,  along  >'ith  the  corresponding 


Figure  3.  Lift  force  phenomenon. 


(a)  Unrestricted  flow  through  tha  bottoa  clearance  at  low  velocltlca 
raaulta  in  downward  lift  force. 

(b)  Unrestricted  flow  through  tha  bottoa  clearance  at  higher  velocltlca 
Increases  the  negatlra  lift. 

(c)  Choking  effect  beglna,  ao  downward  lift  force  decreaaea. 

(d)  Velocltlca  increaae  and  preaaurea  decraaae  at  a faater  rate  over  the 
top  of  the  pipe  aectlon  than  In  the  reatrlcteo  'low  through  the  bottoa 
clearance,  no  the  lift  force  becoaea  poeltive. 

(a)  Upward  lift  forca  lncreaaaa  with  lncreaaing  velocities. 

(f)  Poeltive  lift  reaches  a aazlaua. 

(g)  Poeltive  lift  force  decreases  and  choking  effect  dlalnlahea  with 
decreasing  velocities  of  the  passing  wave  creat. 

(h)  Lift  force  again  becoaea  negative  as  the  flow  through  the  bottoe 
clearance  becoaea  leas  restricted. 

(1)  Unrestricted  flow  through  bottoa  clearance  at  low  velocities  results 
In  downward  lift  force. 

(j)  Lift  force  cycle  la  repeated  aa  the  flow  reverses  with  the  approaching 
wave  ttougit. 


The  traditional  lift  force  equation,  derived  for  unidirectional 
steady- flow  situations,  is  expressed  as  = 1/2  p A u2,  where  is 
the  coefficient  of  lift.  This  equation  has  been  applied  to  wave-induced 
lift  forces,  using  the  horizontal  component  of  the  oscillating  water 
particle  velocity,  u,  in  the  relationship.  The  lift  force  expressed  in 
this  way  assumes  that  the  force  acts  in  one  direction  only  (either  upward 
or  downward)  throughout  the  entire  wave  cycle. 

A pipeline  located  on  the  ocean  floor  with  no  clearance  will 
experience  an  upward  lift  force  throughout  the  entire  wave  cycle, 
increasing  with  the  horizontal  velocities  to  reach  maximum  values  under 
the  crests  and  troughs  of  the  passing  waves,  and  diminishing  to  zero 
as  the  horizontal  velocities  go  to  zero  at  the  point  of  flow  reversal. 

Th: c phenomena  is  described  adequately  by  the  above  lift  force  equation 
with  a positive  coefficient  of  lift  C^. 

A pipeline  located  at  a large  enough  clearance  above  the  bottom  so 
tnat  the  choking  effect  does  not  occur  will  experience  a downward  lift 
force  throughout  the  wave  cycle,  since  the  flow  is  always  faster  through 
the  bottom  constriction  than  over  the  top  of  the  pipeline.  Again,  this 
negative  lift  force  increases  with  the  horizontal  water  particle  veloci- 
ties, reaching  maximum  magnitudes  under  the  crests  and  troughs  of  the 
passing  waves,  and  decreasing  to  zero  as  the  flow  reverses.  This  phe- 
nomenon is  also  suitably  expressed  by  the  traditional  lift  force  equation 
but  using  a negative  coefficient  of  lift. 

These  two  situati  represent  the  extreme  cases  bounding  the  lift 
force  phenomena.  However,  the  choking  phenomenon  will  occur  at  any 
clearance  between  these  two  limits,  and  the  traditional  lift  force 
equation  cannot  be  used  to  accurately  describe  the  forces  exerted  on 
a pipeline.  This  equation  must  be  replaced  by  a model  developed  speci- 
fically for  wave-induced  lift  forces.  The  experimental  results  of  this 
investigation  demonstrate  that  the  largest  wave-induced  lift  forces 
occur  at  these  intermediate  clearances,  where  the  choking  phenomenon 
does  develop. 

Since  the  lift  force  phenomenon  is  repeated  twice  per  wave  cycle 
with  the  reversal  of  the  horizontal  flow  pattern,  the  lift  force  can 
be  described  mathematically  by  a sinusoidal  function  of  twice  the  fre- 
quency of  the  waves.  Jn  addition,  the  mathematical  expression  must 
allow  for  description  of  the  following  lift  force  properties: 
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Figure  4.  Change  in  lift  force  record  for  increasing  bottom  clearance. 
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(a)  The  lift  force  may  be  positive  during  part  of  the  wave  cycle 
and  negative  for  the  rest  of  the  cycle.  The  proportion  of  positive  lift 
to  negative  lift  may  range  from  all  positive  lift  to  all  negative  lift. 

(b)  The  positions  of  the  maximum  values  of  both  the  upward  and 
downward  lift  forces  will  shift  with  respect  to  the  position  of  the 
wave  cycle  as  the  bottom  clearance  is  increased  (for  a given  pipeline 
and  wave  condition) . 

(c)  As  the  clearance  is  increased,  the  maximum  value  of  the  upward 
lift  force  will  decrease,  while  correspondingly  the  maximum  value  of 
the  downward  lift  force  will  increase. 


(d)  When  the  bottom  clearance  is  increased  to  a point  at  which  the 
lift  effect  is  downward  throughout  the  entire  wave  cycle,  further 
increases  in  clearance  will  result  in  decreases  in  the  maximum  magnitude 
of  the  downward  liit  force,  but  without  a shift  in  the  position  of  the 
maximum  lift  force  with  respect  to  the  position  of  the  wave  cycle  over 
the  pipeline . 

A lift  force  equation  of  the  form, 

fL  = 1/2  P A umax2[cos2  (8  - <10  - k],  (4) 

allows  an  adequate  mathematical  description  of  all  the  above  properties 
of  the  wave-induced  lift  force  phenomena.  This  equation  fits  the  experi- 
mental data  reasonably  well  over  the  wide  range  of  conditions  tested. 

The  parameters  involved  in  this  modified  form  of  the  traditional 
lift  force  equation  are: 


Cl  = coefficient  of  lift 

p = mass  density  of  fluid 

A = projected  area  of  pipe  section 

uinax  = maximum  value  of  horizontal  component  of  water 
pai'ticle  velocity  at  center  of  pipe  section  if 
pipeline  was  absent 
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0 = -y-  = position  of  wave  cycle  over  center  of  pipe  section 
with  respect  to  time,  where  T is  the  wave  period 
and  t is  the  time  since  the  last  crest  passed  over 
the  center  of  the  pipe  section  (see  definition 
sketch  in  Fig.  5).  The  wave  crest  corresponds  to 
0=0°  (0  radians)  or  2tt/T  - 0 radians.  The  wave 
trough  corresponds  to  180°  (it  radians)  or 
2irt/T  = it  radians 


<p  = phase  shift  of  maximum  lift  forces  with  respect 

to  wave  cycle 

k * negative  fraction  of  lift  force  cycle 

The  parameter,  k,  represents  the  increase  in  the  magnitude  and 
duration  of  the  negative  lift  forces  acting  on  a pipeline  with  increas- 
ing bottom  clearance,  and  the  corresponding  decrease  in  the  magnitude 
and  duration  of  the  positive  lift  forces.  The  value  of  k varies  from  a 
minimum  of  0 to  a maximum  valxie  of  1.  k = 0 corresponds  to  the  case  of 
a pipeline  lying  on  the  bottom  with  no  clearance,  in  •'nich  the  lift 
forces  are  positive  throughout  the  wave  cycle,  k increases  with 
increasing  bottom  clearance  to  a maximum  value  of  1,  which  corresponds 
to  the  case  of  a pipeline  located  at  a sufficient  clearance  from  the 
bottom  so  that  the  choking  phenomenon  does  not  occur,  and  in  which  the 
lift  forces  are  therefore  negative  throughout  the  wave  cycle. 

The  phase  shift  parameter,  <f>,  represents  the  shift  in  the  position 
of  the  maximum  values  of  both  the  positive  and  negative  lift  forces 
with  respect  to  the  wave  cycle  as  the  bottom  clearance  increases.  The 
value  of  <f>  may  range  from  0°  to  a maximum  value  of  90°.  <f>  = 0°  corre- 

sponds to  the  case  of  a pipeline  located  on  the  ocean  floor  with  no 
bottom  clearance,  in  which  the  lift  forces  are  positive  throughout  the 
wave  cycle  with  maximum  forces  occurring  under  the  crests  and  troughs 
of  the  passing  waves.  <j)  increases  with  increasing  bottom  clearance  to 
a maximum  value  of  90s,  corresponding  to  a pipeline  located  above  the 
bottom  at  a sufficient  clearance  so  that  the  choking  effect  does  not 
occur;  the  lift  forces  are  negative  throughout  the  w rve  cycle  with  maxi- 
mums  occurring  under  the  crests  and  troughs  of  the  waves.  As  defined, 

<j>  = 0*  when  k * 0,  and  = 90°  when  k = 1,  or  vice  versa. 

The  coefficient  of  lift,  Cp,  in  this  form  of  the  lift  force  equation 
will  always  have  a positive  value,  since  negative  values  of  the  lift 
force  are  accounted  for  by  the  value  of  the  parameter,  k.  The  lift 
force  equation  is  shown  graphically  in  Figure  6. 

To  apply  the  lift  force  equation  to  a practical  design  situation, 
values  of  Cp,  k,  and  <j>  must  be  determined  for  a given  set  of  pipeline 
and  wave  conditions  corresponding  to  the  particular  case  under  considera- 
tion. Selection  of  tne  appropriate  values  requires  quantitative  knowl- 
edge of  the  functional  relationships  between  these  parameters  and  the 
wave  conditions,  bottom  clearance,  and  pipeline  size  and  configuration. 

The  development  of  these  relationships  was  the  purpose  of  the  experimen- 
tal part  of  this  investigation. 

In  a real  situation,  a pipeline  on  the  ocean  floor  is  often  laid  over 
an  irregular  bottom,  supported  by  the  high  points  in  the  bottom  topography 
hut  probably  spanning  the  depressed  areas.  In  this  case,  the  pipeline 
must  be  broken  into  component  sections  of  the  same  approximate  bottom 
clearance  for  a separate  analysis  of  each  section.  The  results  of  the 
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analysis  will  yield  the  lift  force  record  (both  magnitudes  and  time 
history)  of  each  separate  component  pipe  section,  which  may  then  be 
integrated  in  the  appropriate  manner  to  determine  the  maximum  wave- 
induced  stresses  exerted  on  the  pipeline  at  any  critical  section. 

This  is  important  because  the  maximum  lift  forces  may  act  upward 
on  a bottom-supported  section  of  a pipeline,  while  acting  downward  on 
the  adjacent  sections  of  the  pipeline  spanning  the  bottom  at  a small 
clearance.  Maximum  values  of  both  the  positive  and  negative  lift 
forces  acting  in  opposite  directions  could  easily  occur  at  the  same 
point  in  the  wave  cycle  (under  the  crests  and  troughs) , thus  exerting 
stresses  on  the  pipeline  twice  as  high  as  would  be  calculated  consid- 
ering any  pipe  section  alone,  or  in  using  some  average  clearance  for 
a long  section  of  the  pipeline. 

4 . Extension  of  Model  to  Higher  Order  Theories. 

The  lift  force  model  (eq.  4)  is  based  on  linear  theory,  assuming 
the  lift  force  phenomenon  is  identical  as  either  the  wave  crest  or 
trough  passes  over  the  pipeline.  Such  a symmetrical  expression  is  not 
flexible  enough  to  consider  slightly  different  kinematics  under  the 
v^ve  crests  and  troughs,  which  are  expressed  in  higher  order  theories. 
Tht'.'.-e  different  kinematics  would,  in  reality,  produce  slightly  different 
lift  forces  under  the  crests  and  troughs  of  nonlinear  waves. 

The  lift  force  model  described  above  was  derived  as  a modification 
of  the  traditional  lift  force  equation  using  linear  wave  theory  to 
express  the  horizontal  water  particle  velocities.  Using  linear  wave 
theory,  the  traditional  lift  force  equation  can  be  expressed  as: 


FL  =■  1/2  CL  p A ^ cos*  (6). 


(5) 


This  equation  was  modified  to  make  it  a suitable  expression  for  wave- 
induced  lift  forces  by  adding  the  phase  shift  parameter,  <j>,  to  account 
for  maximum  lift  forces  occurring  in  places  other  than  the  crest-  and 
trough  in  the  wave  cycle,  and  by  adding  the  parameter,  k,  to  account 
for  positive  lift  forces  during  part  of  the  wave  cycle  and  negative 
forces  during  the  rest  of  the  cycle.  This  modified  equation  fits  the 
experimental  data  very  well  for  all  conditions  tested  in  this  investi- 
gation. 

The  model  was  developed  after  thorough  inspection  of  the  experimental 
data.  For  a given  pipe  diameter  and  wave  condition,  the  force  record 
followed  a sinusoidal  relationship  of  twice  the  frequency  of  the  waves. 

As  the  clearance  increased,  the  maximum  positive  forces  gradually  dimin- 
ished while  continuously  shifting  to  a maximum  of  90°  from  the  wave  crest 
as  the  forces  went  to  zero  (Fig.  4).  At  the  same  time,  the  maximum 
negative  forces  slowly  grew  from  a minimum  value  of  zero  at  a position 
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90®  from  the  wave  crest  and  increased  while  continuously  shifting  posi- 
tions to  reach  a maximum  negative  value  at  a position  180°  from  the 
wave  crest  (Fig.  7,  a). 

Since  a sinusoidal  function  of  twice  the  frequency  of  the  wave 
(sin  20  or  cos  20)  can  be  expressed  as  cos20,  using  the  appropriate  trigo- 
nometric relationships,  and  since  the  lift  force  is  a function  of  the  hori- 
zontal velocity  squared  (u^x  cos  ®)2,  using  linear  wave  theory,  the  lift 

force  equation  was  expressed  as  Ft  » 1/2  C.  p A u 2 [cos2  (0  - <J>)  - k]. 

^ l max 

However,  it  is  the  symmetrical  properties;  of  this  equation  and  linear 
wave  theory  that  allow  this  expression  to  work  so  well.  When  higher 
order  wave  theories  are  applied  to  this  relationship,  problems  due  to 
nonsymmetry  are  encountered.  This  is  easily  seen  by  graphically  compar- 
ing the  transition  from  positive  to  negative  lift  forces  with  increasing 
bottom  clearance  with  this  lift  model,  using  both  linear  and  higher  order 
theories . 

The  horizontal  component  of  the  water  particle  velocity  for  both 
Stokes'  third-order  waves  and  linear  waves  is  shown  in  Figure  8,  along 
with  the  corresponding  lift  forces  on  a pipeline  for  the  two  extreme 
cases  of:  (a)  a pipeline  on  the  bottom  with  no  clearance,  and  (b)  a 
pipeline  with  a large  enough  bottom  clearance  so  that  the  choking  phe- 
nomenon does  not  occur.  By  gradually  shifting  the  linear  theory  lift 
force  ave  for  case  (a)  (no  bottom  clearance)  to  the  right  90°  from 
the  w crest,  while  simultaneously  lowering  it  so  that  the  forces 
become  negative,  the  lift  force  curve  for  case  (b)  is  obtained  (com- 
pare Figs.  7 and  8).  This  same  transformation  of  the  wave  force  record 
ws  served  with  increasing  bottom  clearance  in  the  experimental  data. 

. ‘.ever,  if  this  procedure  is  repeated  with  the  Stoked  third-order 
lift  rce  record,  the  correct  force  record  for  case  (b)  is  not  obtained 
(compare  Figs.  7 and  8).  In  reality,  rather  than  a mere  shift  of  the 
force  record  downward  and  to  the  right  with  increasing  bottom  clearance, 
a simult  neous  transformation  of  the  shape  of  the  lift  force  record 
would  aJiO  occur  for  highly  nonlinear  waves.  This  gradual  transforma- 
tion c he  shape  occurring  simultaneously  with  the  shift  would  provide 
a continuous  change  in  the  lift  force  record  with  increasing  clearance 
between  the  two  limiting  cases  (a)  and  (b)  (Fig.  8). 

However,  the  lift  force  phenomenon  is  not  a direct  function  of  the 
instantaneous  water  particle  velocity  acting  at  the  center  of  the  pipe 
section  if  the  pipeline  was  absent.  Rather,  it  is  a complicated  function 
of  the  asymmetrical  distorted  flow  pattern  and  accelerating  velocity 
field  acting  on  the  pipeline,  which  in  turn  causes  the  choking  phenomenon 
to  occur,  with  the  resulting  change  in  the  relative  differences  in  the 
flow  velocities  and  corresponding  pressure  distribution  over  the  top  and 
bottom  of  the  pipeline.  Boundary  layer  flow  through  the  bottom  constric- 
tion, the  formation  of  a turbulent  jet  and  associated  eddies,  and  a cyclic 
change  in  the  location  of  the  stagnation  point  with  the  accelerating 
velocity  field  further  conplicate  matters.  In  addition,  the  eddies  and 
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increased  turbulence  generated  by  the  jet  may  be  swept  back  through 
the  bottom  constriction  as  the  flow  pattern  reverses  with  the  passing 
waves . 

Because  of  this,  development  of  an  accurate  mathematical  description 
of  the  lift  force  phenomena  for  nonlinear  waves  that  would  cover  the 
conplete  transformation  of  the  lift  force  record  with  increasing  bottom 
clearance,  and  yet  be  flexible  enough  to  allow  application  of  any  higher 
order  theory,  would  be  a formidable,  if  not  iwpossible,  task.  Since  the 
lift  force  model  developed  for  linear  theory  seems  to  fit  the  experi- 
mental data  reasonably  well,  even  for  waves  that  were  obviously  nonlinear, 
it  should  provide  a useful  tool  for  engineering  calculations,  e.  mi  though 
it  may  not  be  flexible  enough  and  theoretically  correct,  to  allow  the  use 
of  higher  order  wave  theories.  The  value  of  the  maximum  horizontal 
velocity,  u„,ax,  can  be  calculated  under  the  wave  crest  using  any  higher 
order  wave  theory;  this  value  can  then  be  used  in  the  linear  lift  force 
model,  possibly  giving  a tetter  approximation  of  the  lift  forces  induced 
oy  highly  nonlinear  waves. 

II.  EXPERIMENTAL  INVESTIGATION 
1 . Experimental  Equipment. 

Model  experiments  were  performed  in  three  different  wave  tanks.  The 
two-dimensional  tests  were  done  in  a 1-foot-wide  wave  channel  in  the 
Hydraulic  Engineering  Laboratory  (HEL)  at  the  University  of  California, 
Berkeley.  The  three-dimensional  tests  were  started  in  the  8-foot-wide 
Naval  Architecture  (NA)  tow  tank,  and  then  continued  in  the  8-foot-wide 
HEL  wave  tank  where  the  majority  of  the  experiments  were  conducted, 
both  located  at  the  Richmond  Field  Station  of  the  University  of  California, 
The  1-foot  wave  channel  is  100  feet  (30.48  meters)  long;  the  8- foot  HEL 
wave  tank  and  NA  tow  tank  are  180  and  200  feet  (54.86  and  60.96  meters) 
long,respecti vely.  All  tests  were  conducted  at  approximately  the  middle 
of  the  tanks.  A Stillwater  depth  of  2 feet  (60.96  centimeters)  was  used 
in  the  two  dimensional  tests,  and  a 3- foot  (91.44  centimeters)  water 
depth  was  used  in  the  three-dimensional  experiments. 

A flapper-type  generator  is  located  at  one  end  of  each  of  the  HEL 
wave  tanks;  the  NA  tow  tank  has  a piston-type  wave  generator.  The  wave 
period  is  controlled  by  varying  the  speed  of  the  electric  motors  which 
drive  the  wave  generators.  A cam  mechanism  with  a variable  stroke  length 
is  connected  between  the  drive  motor  and  the  flapper,  and  the  wave  height 
is  varied  by  changing  the  stroke  length.  A wave  filter,  consisting  of  a 
series  of  vertical  screens,  was  placed  in  front  of  the  wave  generator  in 
the  1- foot-wide  wave  channel  to  smooth  out  any  irregularities  in  the 
generated  waves  due  to  reflections  from  the  flapper.  A permeable  beach 
was  installed  at  the  opposite  end  of  each  of  the  tanks  to  absorb  the 
wave  energy  and  minimize  the  wave  reflections  from  that  end  of  the  wave 
tank. 
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The  wave- induced  forces  on  the  model  pipe  section  were  measured  by 
a wave  force  meter  designed  and  built  by  Al-Kazily  (1972).  A few  modi- 
fications were  made  to  make  the  instrument  more  suitable  for  this  inves- 
tigation. The  same  transducer  unit  was  used  in  all  of  the  experiments, 
but  fittings  of  different  sizes  were  made  to  accommodate  test  cylinders 
of  various  diameters. 

The  force  transducer  consists  of  a strain  bar  mounted  between  two 
supports.  The  model  pipe  section  is  mounted  to  the  strain  bar  in  such 
a way  that  forces  on  the  pipe  induce  bending  stresses  on  the  strain  bar. 
These  forces  are  measured  by  four  strain  gages  mounted  to  the  strain 
bar  at  sections  of  maximum  strain,  with  two  gages  in  compression  and  the 
other  two  in  tension.  The  strain  gages  are  wi-red  in  a Wheatstone  bridge, 
which  is  connected  to  a carrier  amplifier  which  amplifies  the  output  from 
the  strain  gages.  The  signal  is  then  recorded  on  a strip-chart  recorder. 

The  original  strain  gages  were  Bean-type  BAB-13-125DD-120S,  and  were 
mounted  to  the  steel  strain  bar  with  EPY-1S0  two-part  epoy  , and  then 
coated  with  Dow  Coming  Silastic  RTV  silicon  rubber  for  ’ aterproofing. 
Shortly  after  the  beginning  of  the  three-dimensional  test;,  problems 
were  encountered  in  the  operation  of  the  transducer.  The >e  problems 
were  caused  by  the  deterioration  of  the  original  strain  gage  adhesive 
and  coating,  so  new  strain  gages  were  installed  on  the  transducer  unit. 

The  new  gages  were  Micromeasurement- type  EA-06-125AD-120,  bonded  to  the 
strain  bar  with  Micro measurement  M-Bond  610  two-part  strain  gage  adhesive, 
and  then  coated  with  Micromeasuremervt  M-Coat  D and  M-Coat  G for  water- 
proofing protection.  About  halfway  through  the  three-dimensional  tests, 
further  problems  were  encountered  in  the  operation  of  the  transducer  unit, 
probably  due  to  water  leakage  into  the  waterproof  coating.  There  was 
also  evidence  of  corrosion  on  the  steel  strain  bar,  so  it  was  decided 
to  build  a new  force  transducer  using  a stainless-steel  strain  bar  to 
minimize  corrosion,  and  encapsulated  strain  gages  to  minimize  problems 
with  water  leakage.  The  new  strain  gages  were  Micromeasurement-type 
CEA-06-125UW-120.  The  same  strain  gage  adhesive  and  waterproof  coatings 
were  used,  with  Micromeasurement  M-Coat  B along  the  lead  wires  to  mini- 
mize the  change  of  water  "wicking"  along  the  lead  wires  to  the  inside 
of  the  coating  materials. 

The  transducer  mounting  arrangement  was  different  for  the  two- 
dimensional  and  three-dimensional  experiments.  The  test  cylinder  and 
transducer  unit  for  the  two-dimensional  tests  were  mounted  between  two 
support  brackets  on  each  side  of  the  1-foot  wave  channel.  For  the 
three-dimensional  experiments,  the  test  cylinder  3nd  transducer  unit 
were  mounted  between  two  long  dummy  pipe  sections,  which  were  in  turn 
mounted  to  a steel  base.  The  force  meter  and  mounting  arrangement  is 
shown  in  Figures  9 and  10  for  the  two-dimensional  tests,  and  in  Figures 
11  to  15  for  the  three-dimensional  tests. 

A parallel -wire  resistance-type  wave  gage  was  used  to  record  the  waves 
passing  over  the  model  pipe  section.  The  gige  was  mounted  directly  over 
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Figure  14.  Schematic  of  pipeline  model 
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the  center  of  the  test  section,  so  that  the  wave  records  could  be  cor- 
related directly  with  the  resulting  wave-induced  force  record. 

A Brush  dual-strain  gp*w  amplifier  was  used  in  the  experiments, 
with  one  channel  connected  to  the  wave  gage,  and  the  other  channel 
connected  to  the  force  meter.  The  amplifier  was  connected  to  a Brush 
two-channel  rectilinear  writing  recorder  which  continuously  recorded 
the  waves  and  corresponding  wave-induced  forces  on  the  pipe  section 
(Fig.  16). 

An  electronic  digital  data  acquisition  system  (Paulling  and  Sibul, 
1968)  was  used  in  the  three-dimensional  experiments.  The  digitizer  was 
connected  in  parallel  with  the  strain  gage  amplifier  to  record  simulta- 
neously the  wave  and  corresponding  force  data  on  magnetic  tape,  while 
at  the  same  time  the  data  were  being  recorded  continuously  on  the  strip- 
chart  recorder  (Fig.  17).  The  digitizer  sampled  altematingly  from 
both  the  wave  record  and  force  record  at  a rate  of  100  samples  per 
second,  resulting  in  50  samples  per  second  from  each  of  the  two  chan- 
nels. 

2.  Procedure  for  TWo-Dimensional  Experiments. 

a.  Calibration.  Both  the  wave  gage  and  the  force  transducer  were 

calibrated  before  each  set  of  experimental  runs.  The  wave  gage  was 
calibrated  statically  by  raising  and  lowering  the  gage  in  increments  of 
u.ui  root  cent  meters)  and  recording  the  output.  The  force  meter 

was  also  caliorated  statically  by  hanging  weights  in  increasing  equal 
increments  from  a system  of  pulleys  connected  to  the  force  meter  and 
recording  the  output  on  the  strip  chart.  The  force  transducer  was  cali- 
brated in  both  the  upward  and  downward  directions  by  rearranging  the 
pulley  system  and  repeating  the  above  procedure.  The  calibration  method 
is  shown  in  Figure  18. 

b.  Procedure . After  calibrating  the  force  meter,  the  model  pipe 
section  was  lowered  and  fixed  in  a horizontal  position  at  the  desired 
clearance  above  the  bottom  of  the  wave  channel,  with  the  long  axis  of 
the  test  cylinder  parallel  to  the  approaching  wave  crests.  A sliding 
point  gage  was  mounted  to  the  wave  channel  above  the  pipe  section  and 
was  used  to  accurately  set  the  model  pipe  to  the  desired  bottom  clear- 
ance and  aline  the  pipe  section  parallel  to  the  wave  crests.  Once  the 
model  was  in  the  correct  position,  the  mounting  brackets  and  support 
struts  were  clamped  to  the  sides  of  the  wave  channel.  The  force  trans- 
ducer was  mounted  in  such  a way  that  it  was  sensitive  only  to  forces 
acting  in  the  vertical  direction. 

After  the  model  pipe  section  was  mounted  in  position,  the  wave  gage 
was  lined  up  directly  over  the  center  of  the  pipe  section  with  a plumb 
bob  and  then  clamped  in  position.  The  wave  gage  was  then  calibrated  as 
described  above.  The  experimental  arrangement  is  shown  in  Figure  19. 
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Figure  16.  Brush  recording  instruments. 
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The  pipe  model  and  wave  gage  were  mounted  in  a glass-walled  part 
of  the  tank  near  the  middle  of  the  wove  channel  to  facilitate  the  vis- 
ual observation  of  the  phenomenon  being  studied.  For  each  bottom 
clearance  tested,  a series  of  runs  was  made  with  waves  generated  at 
19  different  wave  periods,  covering  a range  of  0.95  to  2.5  seconds. 

Seven  wave  heights  were  generated  for  each  wave  period,  ranging  up  to 
0.34  foot  (10.4  centimeters). 

After  these  runs  were  completed,  the  pipeline  was  set  at  another 
bottom  clearance,  and  the  procedure  was  repeated.  Seven  bottom  clear- 
ances were  tested  for  each  wave  condition,  ranging  from  0.001  foot, 

1/16,  1/8,  3/16,  1/4,  1,  and  2 inches  (0.305,  1.59,  3.18,  4.76,  and 
6.35  millimeters,  2.54  and  5.08  centimeters),  respectively.  The  mini- 
mum clearance  tested  (0.001  foot)  was  that  which  placed  the  pipe  sec- 
tion as  close  to  the  bottom  as  possible  without  touching  the  bottom 
when  the  waves  passed  over  it.  This  was  necessary  to  measure  any  down- 
ward forces  exerted  on  the  pipe  section  due  to  the  wave  action.  The 
2-inch  bottom  clearance  placed  the  pipe  section  far  enough  from  the 
bottom  so  that  the  vertical  lift  forces  were  insignificant. 

These  experiments  were  carried  out  with  a 4-inch-diameter  (10.1b 
centimeters)  test  cylinder.  The  experiments  were  repeated  with  pipe 
sections  of  2-,  2-1j-,  and  3-inch  (5.08,  6.35,  and  7.62  centimeters) 
diameters,  but  only  three  bottom  clearances  were  tested — 0.001  foot, 

1/8  inch,  and  1/4  inch.  The  wave  conditions  covered  the  same  range  of 
wave  heights  and  periods,  but  were  not  quite  as  extensive  in  number. 

In  addition  to  the  vertical  force  measurements,  a series  of  experi- 
ments was  performed  to  measure  the  horizontal  forces  acting  on  the  pipe 
section,  so  that  the  resuJta  t.  wave-induced  force  could  be  determined 
throughout  the  entire  wave  cycle  for  several  of  the  experimental  condi- 
tions tested.  Only  the  4-inch-diameter  test  cylinder  was  used  in  these 
experiments,  since  the  corresponding  vertical  experiments  were  the  most 
extensive  for  the  4-inch  cylinder.  The  horizontal  forces  were  measured 
by  rotating  the  force  transducer  90°  so  that  it  was  sensitive  only  to 
forces  acting  in  the  horizontal  direction.  The  calibration  procedure 
was  the  same  as  described  above  for  the  vertical  force  measurements 
except  that  the  system  of  pulleys  was  rearranged  so  that  the  calibration 
weights  exerted  forces  in  the  horizontal  direction  only. 

All  seven  of  the  bottom  clearances  used  in  the  vertical  experiments 
were  also  used  in  the  horizontal  tests.  The  wave  periods  covered  the 
same  range  as  the  vertical  experiments,  but  only  6 of  the  19  wave  peri- 
ods were  used--0.95,  1.25,  1.5,  1.85,  2.25,  and  2.55  seconds.  Two  of 
the  seven  wave  heights  corresponding  to  each  wave  period  ip  the  vertical 
experiments  were  used  in  the  horizontal  tests. 

The  Stillwater  depth  was  held  constant  at  a depth  of  2 feet  through- 
out the  two-dimensional  tests. 
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3.  Procedure  for  Three-Dimensional  Experiments. 


a.  Calibration.  The  wave  gage  and  force  meter  were  calibrated 
before  each  set  of  experimental  runs.  The  wave  gage  was  calibrated  in 
the  same  manner  as  the  two-dimensional  tests,  but  0.1-foot  (3.05  centi- 
meters) increments  were  used  rather  than  0.05-foot  increments,  since 
larger  waves  were  used  in  these  experiments. 

The  force  transducer  was  calibrated  in  the  upward  direction  in  the 
same  manner  as  the  two-dimensional  tests,  by  hanging  weights  over  a 
pulley  to  a string  attached  to  the  pipe  test  section.  However,  because 
the  three-dimensional  model  was  mounted  to  a base  with  a smal ’ bottom 
clearance,  it  was  impossible  to  calibrate  the  transducer  in  the  down- 
ward direction  by  using  a system  of  pulleys,  since  there  was  no  room 
for  a pulley  between  the  pipe  section  and  the  base  to  which  it  was 
mounted.  Rather,  the  force  meter  was  calibrated  in  the  downward  direc- 
tion by  placing  the  weights  directly  on  top  of  the  center  of  the  sub- 
merged test  section  and  using  the  submerged  weight  of  the  weights  in 
calculating  the  calibration  curve.  Weight  increments  of  50  grams  were 
used  in  calibrating  the  transducer.  The  calibration  method  is  shown  in 
Figure  20. 

b.  Procedure . An  overhead  crane  was  used  to  lower  the  pipeline 
model  and  base  into  the  wave  tank.  The  assembly  was  first  submerged  to 
a depth  of  about  feet  (45.7  centimeters).  The  model  was  tilted  at 
both  ends  to  remove  all  air  bubbles  from  the  system,  and  the  ends  of 
the  dummy  pipe  sections  were  stoppered  to  prevent  waterflow  through  the 
pipeline  model.  The  bottom  clearance  between  the  base  and  the  pipe 
model  was  adjusted  by  placing  spacers  on  the  support  rods  between  the 
base  and  the  dummy  pipe  sections,  and  then  tightening  the  nuts  on  the 
support  rods  above  the  dummy  pipe  sections.  The  test  section  was  tnen 
centered  and  adjusted  carefully  to  the  exact  bottom  clearance  desired 
with  the  aid  of  10  adjusting  screws.  The  calibration  string  was 
attached  to  the  test  section,  and  the  assembly  was  lowered  to  the 
bottom  of  the  tank. 

The  calibration  string  and  pulley  system  was  alined  di : ectly  over 
the  center  of  the  test  section  with  a plumb  bob,  and  the  pul  ey  sup- 
ports were  then  clamped  to  the  sides  of  the  wave  tank.  The  transducer 
was  first  calibrated  in  the  upward  direction,  after  which  the  calibra- 
tion string  was  removed,  and  the  transducer  "was  calibrated  in  the  down- 
ward direction,  as  described  above.  j 

-* 

The  pipeline  model  wa?  positioned  at  the  desired  angle  of  orienta- 
tion on  the  tank  bottom  by  lining  up  one  of  the  long  edges  of  the  model 
base  parallel  to  the  correct  line  marked  on  the  bottom  of  the  wave  tank. 
Lines  were  marked  or  the  tank  bottom  in  15°  increments  from  0°  to  75°, 
where  0°  corresponds  to  a pipeline  parallel  to  the  approaching  wave 
crests.  After  the  model  >vas  calibrated  and  placed  in  position,  the 
wave  gage  was  lined  up  directly  over  the  center  of  the  test  section  with 


a plumb  bob,  clamped  position,  and  then  calibrated  as  described 
above.  The  experimental  arrangement  is  shown  in  Figure  21. 

For  each  bottom  clearance,  six  angles  of  orientation  (0°,  IS®,  30°, 
45®,  60®,  and  75°)  were  tested.  Fifteen  runs  with  different  wave  con- 
ditions were  made  for  each  bottom  clearance  and  orientation  angle. 

These  runs  covered  four  wave  periods  ranging  from  1.4  to  2.6  seconds, 
with  waves  generated  at  four  heights  for  each  period,  ranging  to  a maxi- 
mum of  about  0.7  foot  (21.3  centimeters).  Eight  bottom  clearances  were 
tested,  ranging  from  0.001  foot,  1/16  inch,  1/8  inch,  3/16  inch,  1/4 
inch,  1/2  inch,  1 inch,  and  2 inches. 

The  above  experiments  were  done  using  a 3-inch-diameter  pipeline 
model.  The  tests  were  then  repeated  using  a 2-  and  4-inch-diameter 
pipeline.  The  1-  ar.d  2-inch  clearances  were  not  tested  because  the  lift 
forces  at  these  clearances  proved  insignificant  in  the  previous  tests. 
Also,  the  tests  at  an  orientation  angle  of  75°  were  eliminated,  since 
the  previous  experiments  demonstrated  that  the  vertical  forces  measured 
at  this  angle  were  insignificant,  and  too  small  to  be  measured  with  any 
accuracy.  Aside  from  these  changes,  the  4-inch-diameter  pipeline  was 
tested  at  the  same  bottom  clearances,  orientation  angles,  and  wave  con- 
ditions as  the  3-inch-diameter  model.  The  2-i,ich-diameter  model  was 
tested  at  the  same  bottom  clearances  and  wave  conditions,  but  only  three 
of  the  five  orientation  angles  (0®,  30°,  and  60®)  were  tested. 

The  Stillwater  depth  in  the  wave  tank  was  held  constant  at  3 feet 
throughout  the  three-dimensional  experiments,  but  since  the  base  of  the 
pipeline  model  was  located  2-7/16  inches  (6.19  centimeters)  above  the 
tank  bottom,  the  effective  Stillwater  depth  over  the  pipeline  base  was 
2.797  feet  (85.25  centimeters).  The  definition  sketch  for  the  three- 
dimensional  experiments  is  shown  in  Figure  22. 

4 . Data  Reduction. 

The  wave  force  data  were  taken  on  a two- channel  strip-chart  recorder 
with  the  paper  advancing  at  a speed  of  25  centimeters  per  second.  One 
channel  recorded  the  forces  while  the  other  channel  simultaneous1 y re- 
corded the  wave  surface  profile  directly  over  the  center  of  the  pipeline 
test  section,  thus  allowing  direct  correlation  of  the  two  records. 

The  two-dimensional  experimental  data  were  digitized  manually  using 
a Gerber  digital  data  reduction  system  connected  with  a card  punch  to 
automatically  punch  the  digitized  values  on  computer  cards.  Using  a 
variable  linear  scale,  each  force  record  was  first  ^..  vided  into  20 
equally  spaced  intervals  per  wave,  each  interval  representing  a time 
interval  of  T/20,  where  T is  the  wave  period.  Each  force  record  was 
digitized  at  these  points  over  an  interval  of  two  consecutive  waves 
(beginning  at  the  wave  crest),  thus  giving  -u  > ilues  for  the  analysis 
and  averaging  the  wave  forces  over  two  wave  cy< les. 
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The  points  in  the  force  records  corresponding  to  the  wave  crests 
were  chosen  as  the  origin  (and  end)  of  the  digitized  records.  These 
points  were  determined  by  averaging  the  midpoints  of  three  or  four  hori- 
zontal lines  drawn  through  the  crests  of  the  wave  record  at  several 
elevations  above  the  Stillwater  level  (SWL) . These  midpoints  were 
approximately  identical  except  for  some  of  the  larger,  longer  waves  in 
which  the  peak  of  the  wave  crest  did  not  exactly  coincide  with  the  mid- 
point of  the  zero  crossings  of  the  wave  crest.  A sample  data  record  is 
given  in  Figure  23. 

The  three-dimensional  experimental  da"a  were  handled  differently; 
the  data  were  recorded  on  magnetic  tape  with  an  electronic  digital  data 
acquisition  system.  This  instrument  sampled  alternatingly  from  the  two 
channels  (wave  and  force)  at  a rate  of  100  samples  per  second,  result- 
ing in  50  samples  per  second  from  each  channel. 

The  origin  at  the  wave  crest  and  the  v,.ve  period  were  determined 
from  the  digitized  wave  records,  rathe’*  than  directly  from  the  strip- 
chart  records.  Since  positive  readings  of  che  wave  profile  corres- 
ponded to  the  crest  ard  negative  readings  corresponded  to  the  trough, 
the  point  of  origin  of  the  wave  crest  was  determined  by  taking  the  mid- 
point of  the  positive  readings  between  zero  crossings  on  the  wave 
profile.  The  crest  was  thus  defined  as  the  data  point  closest  to  the 
midpoint  of  the  zero  crossings.  The  wave  period  was  determined  from 
the  number  of  readings  between  two  successive  crests,  since  there  was 
a time  interval  of  1/50  second  between  each  reading.  Thus,  the  wave 
period  was  determined  tn  the  nearest  0.02  second. 

The  origin  of  the  force  record  wac  taken  as  the  force  reading  cor- 
responding to  the  defined  origin  at  the  center  of  the  wave  crest  surface 
profile,  in  reality,  there  was  a small  time  lag  of  1/100  second  between 
the  wave  profile  readings  and  the  corresponding  force  readings.  This 
small  timelag  was  ignored  in  the  analysis,  since  it  was  felt  that  the 
accuracy  of  the  defined  origin  at  the  wave  crest  was  only  good  to  the 
nearest  1/50  second,  the  time  interval  between  successive  readings  of 
the  wave  record. 


Only  one  wave  cycle  was  used  for  analysis  of  the  electrons cally 
digitized  data.  Since  the  data  were  on  magnetic  tape,  it  was  impossible 
to  determine  that  two  successive  waves  had  exactly  the  same  period  and 
height  until  after  the  calculations  were  completed  on  the  computer. 

Thus,  if  the  waves  had  slightly  different  periods,  the  time  phase  cor- 
relation of  the  corresponding  force  readings  would  be  slightly  in  error 
when  taken  over  two  wave  cycles.  In  addition,  since  the  accuracy, 
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more  readings  per  wave  cycle  than  the  manual  digitizing  method,  a suf- 
ficiently large  number  of  force  readings  could  be  obtained  in  one  wave 
cycle . 
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An  estimation  of  the  accuracy  of  the  experimental  measurements,  along 
with  the  sources  of  error,  is  tabulated  in  the  Table. 

A least  squares  analysis  was  performed  on  the  digitized  force  data 
to  calculate  the  parameters,  Cp,  4>.  CM,  and  C[),  of  the  vertical  wave 
force  equation,  and  the  coefficients,  and  Cp,  corresponding  to  the 
horizontal  wave  force  equation.  Using  this  approacn,  values  of  the  wave 
force  parameters  that  best  fit  the  force  data  throughout  the  entire  wave 
cycle  can  be  determined.  These  values  were  then  substituted  back  into 
the  wave  force  equation  to  calculate  the  force  over  a complete  wave  cycle, 
thus  allowing  comparison  of  the  results  with  the  original  data.  The 
least  squares  analysis  is  given  in  Appendix  A.  The  computer  programs 
used  for  the  analysis  are  given  in  Appendixes  B,  C,  and  D;  the  tabulated 
results  of  the  analysis  are  in  Appendixes  E,  F,  and  G.  Examples  of  the 
computer  output  showing  comparison  of  the  original  data  with  the  forces 
calculated  using  the  results  of  the  least  squares  analysis  are  given  in 
Figures  24  and  25. 


III.  RESULTS  AND  DISCUSSION 
1 . Resultant  Force  Through  Wave  Cycle. 

Both  horizontal  and  vertical  force  measurements  were  made  for  some 
test  conditions  in  the  two-dimensional  experiments  using  the  4-inch- 
diameter  cylinder.  The  resultant  force  throughout  the  wave  cycle  could 
thus  be  determined  for  these  conditions.  Figures  26  to  32  show  the 
resultant  force  plotted  for  each  bottom  clearance  under  the  same  wave 
condition,  a period  of  1.85  to  1.86  seconds  and  a wave  height  of  0.24 
to  0.25  foot  (7.32  to  7.62  centimeters).  Values  from  the  corresponding 
horizontal  and  vertical  force  records  were  plotted  at  20  evenly  spaced 
intervals  (18°)  through  each  wave  cycle.  The  forces  were  plotted  for 
two  consecutive  wave  cycles  to  indicate  the  degree  of  scatter  in  the 
data.  A rectangle  was  drawn  at  each  plotted  point  to  illustrate  the 
horizontal  and  vertical  range  of  the  force  data  over  the  two  wave  cycles, 
and  an  envelope  curve  was  drawn  over  these  points. 

Examination  of  these  plots  as  a group  (Fig.  33)  shows  the  transition 
of  the  resultant  wave-induced  force  with  increasing  clearance  for  the 
given  wave  condition  (T  = 1.85  to  1.86  seconds,  11  * 0.24  to  0.25  foot). 
The  vertical  component  of  the  wave  force  is  dominated  by  the  lift  force, 
while  the  horizontal  component  of  the  resultant  force  is  due  to  the  iner- 
tial and  drag  forces,  with  the  inertial  forces  predominating  for  the 
experimental  conditions  tested. 

For  the  smallest  clearance  (0.001  foot),  in  which  the  pipeline  is 
almost  in  contact  with  che  bottom,  the  resultant  force  attains  a maximum 
upward  value  under  the  crests  and  troughs  of  the  passing  waves.  The 
total  wave  force  acts  in  the  upward  (positive)  direction  throughout  the 
complete  wave  cycle,  except  for  small  downward  forces  in  the  vicinity  of 
90°  and  270°,  where  tnt*  horizontal  flow  reverses. 
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Wave  force  5 percent,  except  for  data  taken  at  Stability  of  amplifier  with  respect  to 

largest  orientation  angles  (60°  and  75°)  force  calibration 
in  the  three-dimensional  tests,  which 
are  accurate  to  within  10  percent 


e«LQ*tA  0V4*/h  >1  tM*lt«|W4VH  •■••••CIV  » 

VitfCI  MtITTI  •AveklFfl  CktA*l*f*  *"£<**  A I 

i**i  ..*»  «*N  i.h»  .m  .it*  .»»> 


TOtIL  »W  »Q 

.fSJK. 

ihi* 

I.1NU4 

COM 

I.MMM 

» IMA/ AIM/ 

• IMAAA 

*4 

A* 

VI 

iCIMM 

.IMtll 

-•IMAM 

-••At AAA 

• 1 NMl 

AM 

* a# 

• IMM* 

•CAAAM 

•5?i 

ftB. 

CM 

A*  Alt 

wfiu 

,iiu  «!?}• 

TAAMLAt 

.Mill 

* 

+ 

• 

0 

1 
+ 

-.tun 

• A 

• 

A 

••lull 

-•MM* 

AA 

• | 

A 

••It** 

••aMTt 

•• 

• 1 

i 

•.(MM 

•t(KM 

-.M»l* 

• A 

• 1 A 

••(MN 

••Mall 

•••Mil 

A A 

• 1 1 

*• 

• 1 A 

-.•MM 

•.MIN 

•MM* 

• 

• A 

f.IMM 

•«I*H 

•A 

I.  A 

-♦f**A* 

«.*mi 

• MM* 

A A 

At  • 

-.'UW> 

•DIM 

A 

-.MAM 

M 

A 

•>hmv 

■•MM} 

• 

A 

•»S«)M 

-.(»}) 

*♦ 

• 1 

•ilttH 

••IMM 

- .At  1 A* 

A* 

• 1 

••urn 

-•MM* 

• 

♦ 

*1 

•.Hill 

-•JIM* 

-•MM* 

• 

• 1 

•HIM 

-.MIM 

•MM* 

A 

••WIT 

••MM 

•urn 

A 

• A# 

• IUN 

•*••** 

• MM* 

A 

• A 

•HIM 

•fMW 

A 

• ♦ 

•HIM 

•IAIN 

• MCTC 

A 

. • 

• IMM 

•WM 

-.M*n 

A 

• t 

•I 

• *•**« 

-•MAM 

A 

• 1 

• 

*MMI 

dMtl 

-.Mm 

A 

• 1 ♦ 

asm 

• MW* 

-•M*«* 

A 

• 1 • 

'•!*** 

■HIM 

•NMl 

A 

• 

A 

• 

-••MM 

•»JHf 

♦♦  A 

••*•?*« 

-•MCA* 

♦ 

• IA 

••Iltfl 

-•Ail** 

*♦ 

• 1 A 

'•IMM 

-•IMM 

-•MM* 

• 

• 1 A 

•nVlM 

«.fW* 

♦ 

• 1 A 

-<um 

• •**•*> 

• M «T 

M 

• A 

•»»TMT 

•mi« 

*• 

1.  A 

-.mi* 

•.9MM 

•MSI* 

• 

• 

t 

-.IMII 

•*tM4 

.tmi 

•A 

• 

A 

* ■ measured  force 

4*  - calculated  force 

* - measured  force  - calculated  force 
$ ■ wave  surface 


J . 


’ISV 


?.u: 


ItltSlH 

ftP 

-iftftftftftl , 

•i4!Si 


-“cr1 

•< 

• iTPPPft 

.JSU 


COftft 

t<  m 


* ■ Beasured  force 

+ * calculated  force 

• » measured  force  - calculated  force 


Figure  25.  Example  of  computer  output  for  horizontal  lea^t  squares  analysis. 
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As  tht  clearance  is  increased  t.o  \/16  inch,  the  maximum  upward 
forces  decrease  in  magnitude,  and  also  occur  slightly  later  in  the  wave 
cycle.  At  the  same  time,  the  downward  forces  increase,  reaching  tneir 
maximum  values  at  approximate] > 90°  and  270s. 

Further  increases  in  the  bottom  clearance  produce  a continuous  shift 
of  the  positions  of  both  the  maximum  upward  and  maximum  downward  forces 
later  in  the  wave  cycle.  Simultaneously,  the  forces  become  downward 
rather  than  upward  for  a larger  part  of  the  cycle.  At  the  same  time, 
the  vertical  components  of  the  wave  force  under  the  crests  and  troughs 
become  negative  and  increase  in  the  downward  direction,  while  the  negative 
forces  at  90°  and  270°  gradually  decrease  to  zero. 

At  a 1-inch  clearance,  the  resultant  force  acts  downward  throughout 
almost  the  complete  wave  cycle,  with  maximum  downward  forces  occurring 
under  the  crests  and  troughs  of  the  passing  waves.  The  vertical  forces 
are  zero  at  90°  and  770°,  the  positions  of  the  maximum  horizontal  iner- 
tial forces.  However,  the  lift  effect  is  not  very  large  for  the  1-inch 
clearance.  The  resultant  force  plot  for  the  2-inch  clearance  shows  that 
the  lift  effect  is  still  present,  but  is  relatively  small,  even  in  com- 
parison to  the  small  vertical  inertial  forces. 

At  a s]ighrl>  larger  clearance,  the  lift  effect  will  disappear,  and 
the  vertical  forces  will  be  due  almost  entire  1)  to  the  inertial  forces, 
since  the  vertical  drag  forces  are  negligible  near  the  bottom.  At  this 
clea-cance,  the  inertial  force  will  act  upward  under  the  trough  and  down- 
ward under  the  crest,  so  the  resultant  force  plot  will  take  the  form  of 
an  approximately  symmetrical  ellipse.  This  condition  is  shown  in 
Figure  34  for  a smaller  wave  period  (1.23  seconds),  with  a 1-inch  bottom 
clearance.  The  ellipse  is  distorted  slightly,  due  to  the  smalx  drag 
forces  acting  in  the  horizontal  direction,  90°  out  of  phase  with  the 
larger  horizontal  inertial  forces. 

The  horizontal  components  of  the  resultant  wave  force  are  also 
affected  by  the  proximity  of  the  bottom  bo>mdary.  Although  the  horizon- 
tal water  particle  velocities  and  acceleration::  increase  with  distance 
above  the  bottom,  the  corresponding  horizontal  drag  ind  inertial  forces 
art  larger  when  the  p.  pe  is  close  to  the  bottom  than  when  it  is  located 
above  at  larger  clearances. 

Figures  35,  36,  and  37  show  the  resultant  force  plots  at  both  large 
and  small  bottom  clearances,  for  a wave  with  a period  of  0.95  to  0.96 
second  and  - height  of  0.24  to  0 25  foot.  Because  the  wave  peiiod  is 
small,  the  horizontal  excursions  of  the  water  particles  at  the  bottom 
and  the  duration  of  the  horizontal  flow  are  tco  small  foi  the  lift 
effect  to  develop.  So  the  forces  acting  in  both  the  horizontal  and  ver- 
tical directions  are  mostly  inertial,  with  a small  drag  component  in  the 
horizontal  direction.  The  resultant  force  plots  therefore  take  the  form 
of  an  ellipse. 


Resultant  fo-ce  through  wave  cycle  for  0.001-foot  clearance 
0.95-second  period,  and  0.24- foot  height. 


0.96-second  period,  and  0.25-foot  height. 


However,  the  horizontal  components  of  the  resultant  forces  are  larger 
at  the  smallest  bottom  clearances,  even  though  the  lift  phenomenon  is 
absent.  The  presence  of  the  bottom  boundary  produces  an  asymmetric  flew 
field  around  the  pipeline.  The  resulting  velocities  and  accelerations  of 
the  water  particles  over  the  pipe  section  are  thus  modified  by  the  presence 
of  the  boundary,  and  the  associated  horizontal  forces  are  larger  than  they 
would  be  if  subject  to  the  same  kinematics  in  the  absence  of  the  boundary. 
'Hie  increased  hor’  aontal  forces  on  pipelines  located  close  to  the  bottom 
are  reflected  in  increased  values  of  the  coefficients  of  mass  and  drag, 
and  Cp. 

2 . Orientation  Angle  Considerations. 

The  coefficient  of  lift  calculated  in  the  least  squares  analysis  of 
the  experimental  data  was  computed  using  two  alternative  approaches 
(Fig.  38):  (a)  the  total  horizontal  water  particle  velocity  in  the 

direction  of  wave  advance,  with  the  projected  area  of  the  pipeline  in  the 
plane  perpendicular  to  the  direction  of  wave  advance;  and  (b)  only  the 
component  of  the  horizontal  water  particle  velocity  perpendicular  to  the 
pipeline  axis,  with  the  projected  area  in  the  plane  parallel  to  the  pipe- 
line axis. 

After  tabulating  the  data  from  the  three-dimensional  experiments,  it 
became  apparent  that  the  second  method  gave  consistent  values  of  the 
coefficient  of  lift  for  all  angles  of  orientation.  In  contrast,  the  values 
of  CL  obtained  using  the  first  method  gave  values  that  were  low,  and  which 
decreased  with  increasing  angles  of  orientation  (where  0°  corresponds  to 
a pipeline  parallel  to  the  wave  crests) . 

Relationships  between  the  coeffici-  at  of  lift.  Cl,  and  the  parameters, 

<J)  and  k,  of  the  lift  force  equation  were  the  same  for  all  angles  of  orien- 
tation when  Cl  was  calculated  considering  only  the  component  of  the  hori- 
zontal velocities  perpendicular  to  the  pipeline  axis. 

In  uddition,  relationships  involving  any  of  the  parameters  of  the  lift 
force  equation  (Cl»  4>»  or  k)  and  various  dimensionless  parameters  defining 
the  wave  and  pipeline  conditions  were  consistent  for  all  angles  of  orien- 
tation whe.-  the  horizontal  water  particle  velocity  acting  on  the  pipe 
section  was  treated  by  considering  only  the  conponent  perpendicular  to  the 
pipeline,  and  completely  ignoring  the  parallel  component. 


Thus,  the  results  of  this  investigation  show  that  the  modified  lift 
force  equations  presented  in  this  report  can  be  applied  to  pipelines 
located  at  any  angle  of  orientation  with  respect  to  the  wave  crests. 

However,  only  the  component  of  the  horizontal  water  particle  velocity 
perpendicular  to  the  pipeline  axis  should  be  considered  as  contributing 
to  the  wave-induced  lift  force  acting  on  the  pipeline.  Using  this  approach, 
the  parameters,  Cl,  <j>,  and  k,  lefinirg  the  lift  forces  exhibit  the  same 
quantitative  relationships  bet  'en  the  various  dimensionless  parameters  de- 


3.  Interrelationships  Between  Cl,  $,  and  k. 

4>  and  k were  defined  as  varying  from  0°  to  90°  and  0 to  1,  respec- 
tively, with  increasing  clearance.  4>  = 0°  and  k * 0 correspond  to  the 
case  of  a pipeline  in  contact  with  the  bottom  (no  clearance),  while  the 
maximum  values  of  <t>  - 90°  and  k = 1 correspond  to  the  case  of  a large 
enough  clearance  so  that  the  choking  phenomenon  does  not  occur  at  any  time 
throughout  the  wave  cycle.  Since  a simultaneous  increase  of  both  parameters 
was  noted  in  the  data  for  increasing  clearance  between  the  two  limiting 
cases,  it  was  suspected  that  a direct  relationship  may  exist  between  4>  and 
k.  Such  a relationship  was  found,  as  shown  in  Figure  39.  The  same  rela- 
tionship held  for  all  three  pipe  diameters  tested,  regardless  of  the  ori- 
entation angle,  indicating  that  the  relationship  was  independent  of  these 
two  factors,  and  was  thus  valid  for  any  pipeline  configuration  in  which 
the  lift  effect  was  present. 

In  this  plot  and  the  ones  that  follow,  the  data  for  orientation  angles 
from  0°  to  30°  were  plotted  for  each  pipe  diameter,  without  differentiating 
the  data  corresponding  to  each  angle.  The  relationships  shown  were  found 
to  be  valid  regardless  of  the  angle  of  orientation,  provided  the  data  were 
handled  as  discussed  above  (using  the  component  of  the  horizontal  velocity 
perpendicular  to  the  pipeline  axis).  The  data  corresponding  to  each  pipe 
diameter  are  distinguished  by  using  different  plot  symbols.  The  same 
relationships  hold  for  orientation  angles  of  45°,  but  these  data  were  not 
plotted  in  order  to  minimize  scatter  so  that  differences  between  the  pipe 
diameters  could  be  detected  more  easily.  In  general,  the  same  relation- 
ships held  for  orientation  angles  up  to  60°.  But  in  some  cases,  the  lift 
effect  was  negligible  at  high  orientation  angles,  so  the  values  of  the 
associated  parameters  (C^,  <j>,  and  k)  were  less  accurate.  Thus,  plotting 
all  of  the  data  corresponding  to  the  larger  orientation  angles  would  intro- 
duce additional  scattert  obscuring  the  valid  relationships  which  were 
consistent  when  the  lift  forces  were  significant. 

A relationship  was  found  between  the  coefficient  of  lift,  C^,  and  the 
parameters,  <j)  and  k (Figs.  40  and  41).  Cl  appears  to  be  better  correlated 
with  k than  with  4>.  Note  that  for  minimum  values  of  k and  <J>,  corresponding 
to  the  case  of  a pipeline  in  contact  with  the  bottom,  the  value  of  Cl  is 
approximately  4.5.  This  value  is  of  interest,  since  it  agrees  with  the 
potential  flow  solution  (Cl  = 4.495)  for  the  value  of  the  coefficient  of 
lift  for  a circular  cylinder  in  contact  with  a plane  wall,  subject  to  an 
inviscid  steady  flow  (Yamamoto,  Nath,  and  Slotta,  1973). 

Maximum  values  of  Cl  occur  at  approximately  k = 1/2,  corresponding  to 
maximum  lift  forces  that  are  equal  in  both  the  upward  nd  downward  direc- 
tions. The  average  value  of  the  coefficient  of  lift  at  this  point  is  about 
9.0,  with  values  extending  up  to  about  10.5.  These  maximum  values  of  Cl 
are  attained  at  approximately  <J>  = 25°  to  30°  in  the  <j>  versus  Cl  plot. 

Since  the  coefficient  of  lift.  Cl,  defines  the  confined  magnitude  of 
both  the  positive  and  negative  lifts,  it  can  be  separated  into  two  parts: 

(a)  the  part  defining  the  magnitude  of  the  positive  lift,  Cl(1-10,  and 
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(b)  the  part  defining  the  magnitude  of  the  negative  lift,  Ci(k)  . The 
quantities,  C^(l.-k)  and  C^Ck) , can  be  referred  to  as  the  effective  posi- 
tive coefficient  of  lift  and  the  effective  negative  coefficient  of  lift, 
respectively.  Since  Cl  = 9.0  for  k = 1/2,  both  C^fl-k)  and  C^(k)  are 
equal  to  4.5  at  this  point.  This  means  that  the  lift  forces  can  reach 
the  same  maximum  magnitude  in  both  the  upward  and  downward  directions  as 
are  attained  in  the  upward  direction  only  for  the  same  pipe  in  contact 
with  the  bottom  (where  C^fl-k)  = 4.5,  but  C^fk)  = 0). 

The  effective  positive  and  negative  coefficients  of  lift  are  plotted 
versus  both  <j>  and  k in  Figures  42  to  45.  Again,  the  correlations  are 
much  better  with  k than  with  <f>.  The  average  value  of  C^Cl-k)  drops  only 
slightly  between  k = 0 and  k = 1/2,  but  for  values  of  k greater  than  1/2, 
the  effective  positive  coefficient  of  lift  drops  rapidly  to  a value  of 
0 when  k = 1 . 

The  average  value  of  C^k)  increases  with  k until  it  reaches  a maximum 
value  of  about  6.0  when  k = 0.75,  and  then  decreases  to  about  4.5  when  k * 1. 
Individual  maximum  values  of  C^k)  attain  values  slightly  greater  than  7.0 
in  the  vicinity  of  k = 0.75.  But  even  the  average  maximum  value  of  6.0  for 
the  effective  negative  coefficient  of  lift  indicates  that  the  downward 
lift  forces  may  attain  maximum  values  33  percent  greater  than  the  maximum 
possible  lift  forces  acting  in  the  upward  direction.  Maximum  values  of 
C,  (k)  corresponds  to  a value  of  <J>  of  about  45°,  which  is  half  way  through 
the  phase  shift  cycle. 

The  potential  flow  theory  gives  a value  of  Cl  = 4.495  for  zero  bottom 
clearance,  with  a discontinuous  jump  to  very  high  negative  values  of  Cl 
for  a very  small  clearance  (Yamamoto,  Nath,  and  Slotta,  1973) . In  the 
potential  flow  solution,  the  value  of  depends  only  on  the  relative 
clearance;  i.e,,  the  ratio  clearance-diameter.  The  coefficient  of  lift  is 
negative  whenever  the  pipe  is  not  in  contact  with  the  bottom,  and  its 
magnitude  decreases  as  the  relative  clearance  is  increased. 

Although  the  potential  flow  solution  appears  to  work  reasonably  well 
when  a pipeline  is  touching  the  bottom,  this  approach  does  not  work  when 
there  is  a small  clearance.  This  is  because  viscous  effects  are  very 
important  for  the  flow  through  the  narrow  bottom  clearance  constriction. 

The  choking  phenomenon  limits  the  maximum  flow  velocities  and  corresponding 
pressure  drops  on  the  bottom  side  of  the  pipeline,  thereby  limiting  the 
maximum  possible  downward  lift  forces. 

The  results  of  this  investigation  indicate  that  the  effective  negative 
coefficient  of  lift,  (^(k^,  can  attain  a maximum  value  of  only  7.0.  This  is 
muen  less  than  the  values  of  Cl  suggested  for  small  relative  clearances  by 
potential  flow  theory.  The  coefficient  of  lift  is  obviously  not  a simple 
function  of  relative  clearance,  since  for  a given  clearance  and  diameter, 
both  the  lift  effect  and  the  coefficient  of  lift  will  vary  with  the  wave- 
induced  flow  conditions.  For  the  smallest  relative  clearances,  the  positive 
lift  forces  were  larger  than  the  negative  lift  forces,  especially  where 
the  horizontal  water  particle  velocities  and  excursions  were  high. 
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Effective  negative  coefficient  of  lift  versus 


Figure  45.  Effective  negative  coefficient  of  lift  versus 


The  largest  negative  lift  forces  do  not  occur  at  clearances  where  the 
choking  effect  is  absent  (corresponding  to  k = 1 and  <f>  = 90°).  Rather, 
the  largest  values  of  the  effective  negative  coefficient  of  lift  correspond 
to  values  of  4>  = 45°  and  k = 0.75.  Interestingly,  when  k = 1 and  = 90° 
where  the  positive  lift  forces  have  decreased  to  zero  and  the  choking 
effect  does  not  develop,  the  i.nxi mum  effective  negative  coefficient  of 
lift  is  approximately  4.5,  the  same  magnitude  as  the  potential  flow  solu- 
tion for  the  positive  coefficient  of  lift  for  zero  bottom  clearance. 
However,  as  the  bottor.  clearance  is  increased  further,  k and  <{>  remain 
at  1 and  90°,  respectively,  while  the  effective  negative  coefficient  of 
lift  decreases  to  zero  (with  the  diminishing  lift  forces). 

The  significance  of  these  results  is  easily  seen  by  following  these 
relationships  for  a given  pipe  and  wave  as  the  pipeline  is  raised  from  the 
bottom,  and  k goes  from  0 to  1.  In  the  interval  from  k = 0 to  1/2,  the 
magnitude  of  the  maximum  upward  lift  forces  remains  the  same,  which  is 
approximately  equal  to  the  potential  flow  solution  for  a cylinder  in 
contact  with  the  bottom  (Cg  = 4.5).  However,  at  the  same  time,  the  nega- 
tive lift  forces  increase  continuously,  reaching  a magnitude  equal  to  the 
positive  lift  forces  at  k = 1/2  (Cg(l-k)  = Cj(k)  = 4.5).  Simultaneously, 
there  is  a shift  in  the  positions  of  both  the  maximum  positive  and  nega- 
tive lift  forces,  since  4>  increases  from  0°  to  30°. 

In  the  interval  k = 1/2  to  1,  the  maximum  positive  lift  rorces 
continuously  decrease  to  zero.  At  the  same  time,  the  maximum  negative 
lift  forces  increase  to  reach  a maximum  value  at  k - 0.75  (where  C^k)  = 

6.  or  7.),  and  then  decrease  back  to  a maximum  corresponding  to  CgOO  = 

4.5  at  k = 1.  The  point  of  maximum  negative  lift  corresponds  to  4>  = 45°, 
the  midpoint  of  the  phase  shift  cycle. 

The  phase  shift  of  the  maximum  lift  forces  is  only  half  as  much  in  the 
interval  k = 0 to  1/2  (where  4>  goes  from  0°  to  30°)  as  in  the  interval 
k = 1/2  to  1 (where  <j>  goes  from  30°  to  90°)  . 

At  k = 1,  only  negative  lift  forces  exist,  and  these  go  to  zero  as 
the  bottom,  clearance  is  increased  further. 

All  of  the  above  interrelationships  between  <j>,  k,  C^,  CgCl-k),  and 
Cg(k)  were  the  same  for  all  pipe  diameters  tested,  regardless  of  the  angle 
of  orientation  (provided  that  was  calculated  considering  only  the  com- 

ponent of  the  horizontal  velocity  perpendicular  to  the  pipeline  axis). 

Thus,  for  the  range  of  conditions  tested,  these  interrelationships  were 
independent  of  the  scale  and  configuration  of  the  pipeline.  Also,  there 
is  no  mention  of  the  wave  conditions,  which  indicates  the  interrelation- 
ships are  independent  of  the  wave  conditions  as  well. 

The  relationships  between  the  parameters,  C^,  4>i  and  k,  defining  the 
lift  force  equation  are  useful,  since  if  either  4>  or  k is  known,  the 
other  two  parameters  can  be  determined.  All  that  is  needed  is  a rela- 
tionship between  <}>  or  k and  the  wave  and  pipeline  conditions. 
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There  appears  to  be  a better  correlation  between  k and  the  parameters 
involving  C,(l-k),  and  C^fk))  than  between  the  analogous  relation- 

ships using  <j>,  so  tlie  former  relationships  should  be  used.  Also,  in  com- 
paring the  plots  of  C^Cl-k)  versus  k (Fig.  42)  and  C^(k)  versus  k (Fig.  44), 
the  scatter  appears  minimal  in  the  plot  with  CL(k)  for  the  interval  of  k 
between  0 and  1/2.  For  the  interval  or  k between  1/2  and  1,  the  scatter 
is  much  less  on  the  plot  between  CL(l-k)  and  k.  Therefore,  it  is  suggested 
that  when  determining  a value  of  for  a given  value  of  k,  the  plot  of 
CL(k)  versus  k be  used  for  values  of  k less  than  1/2  (except  for  k close 
to  0),  and  the  plot  of  CL(l-k)  versus  k be  used  for  values  of  k greater 
than  1/2  (except  for  k close  to  1)  (see  Fig.  46).  For  k close  to  0,  it 
can  be  assumed  that  = 4.5.  However,  for  k ~ 1,  the  value  of  Cl  can 

vary  from  about  4.5  to  zero,  since  as  the  clearance  is  increased  from  the 
point  where  <f>  = 90°  and  k = 1,  both  <+>  and  k remain  at  their  maximum  values 
of  90°  and  1,  respectively,  while  the  lift  effect  diminishes  to  zero. 

When  the  above  relationships  between  <p,  k,  Cl»  CL(l-k),  and  C^(k)  are 
plotted  for  only  the  4-inch-diameter  pipe  model,  the  scatter  is  reduced. 
Although  the  data  for  all  three  diameters  completely  overlap  (showing  the 
same  relationships  hold  for  all  diameters),  the  amount  of  scatter  increases 
with  the  smaller  diameter  models.  This  is  because  the  data  extend  to 
higher  relative  clearances  (clearance-diameter)  for  the  smaller  diameter 
models  than  the  corresponding  data  for  the  4-inch- diameter  model,  since  all 
models  were  tested  at  the  same  actual  clearances. 

Since  the  lift  effect  diminishes  at  high  values  of  the  relative 
clearance,  the  lift  forces  on  the  smaller  diameter  models  at  the  largest 
bottom  clearances  were  very  small  in  many  cases.  This  is  especially  true 
for  the  smaller  waves  and  highe .*  orientation  angles,  where  the  horizontal 
velocities  perpendicular  to  the  pipeline  were  very  low.  In  such  cases, 
the  lift  forces  were  often  insignificant,  so  the  values  of  Cg,  4>»  and  k 
calculated  from  the  least  squares  analysis  were  not  as  accurate. 

In  addition,  as  the  lift  forces  decrease  with  high  relative  clearances, 
eddy-induced  forces  may  approach  the  magnitude  of  the  lift  forces,  thus 
introducing  furtner  error  in  the  calculated  values  of  C^,  <J>,  and  k. 

The  lift  forces  were  generally  significant  for  all  clearances  tested 
using  the  4-inch-diameter  pipe  section,  and  since  the  measured  forces  were 
larger,  the  experimental  error  involved  in  measuring  them  was  less  than 
for  the  smaller  diameter  models. 

Because  of  this,  the  data  taken  for  very  large  bottom  clearances  were 
not  included  in  the  plotted  relationships.  For  higher  clearances,  values 
of  k and  <p  equal  to  1 and  90°,  respectively,  would  be  expected,  since  the 
choking  phenomenon  would  not  occur  throughout  the  vave  cycle.  However, 
as  the  clearance  is  increased,  the  lift  effect  diminishes,  resulting  in 
decreasing  values  of  the  coefficient  of  lift. 
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(k)  versus 


If  such  data  were  included  in  the  plots  of  versus  k and  41,  values 
of  Ct  ranging  from  0 to  the  maximum  values  shownLin  Figures  40  and  41 
would  be  present  in  the  vicinity  of  k * 1 and  4>  » 90®  in  the  respective 
plots.  The  sane  applies  to  the  plots  of  C^(k)  versus  k and  4*. 

These  trends  were  observed  in  the  data  taken  for  the  largest  botton 
clearances  (1  and  2 inches) . However,  since  these  lift  forces  were  so 
small,  a significant  amount  of  error  could  be  introduced  into  the  cal- 
culated values  of  Ci , 4>»  and  k br cause  of  the  presence  of  eddy- induced 
forces,  as  discussed  above.  Therefore,  these  data  were  omitted  from 
the  plotted  relationships,  since  errors  in  4>  or  k corresponding  to  low 
values  of  CL  would  produce  considerable  scatter,  obscuring  the  valid 
relationships  shown. 

4.  Relationships  Between  4>  and  k and  Parameters  Defining  the  Wavs  and 

Pipeline  Conditions. 

To  use  the  above  relationships  between  4>>  and  k to  determine 
the  wave- induced  lift  forces  acting  on  a pipeline,  either  4>  or  k must 
be  known.  Thus,  a value  of  one  of  these  parameters  must  be  determined 
from  relationships  of  4>  or  k with  the  wave  conditions  and  pipeline  con- 
figuration. 

The  lift  force  phenomenon  is  a function  of  the  following  variables: 

(a)  Pipeline  configuration 

(1)  Diameter 

(2)  Clearance 

(3)  Orientation  angle 

(b)  Fluid  properties 

(1)  Density 

(2)  Viscosity 

(c)  Wave- induced  flow  conditions 

(1)  Maximum  horizontal  water  particle  velocity  perpendicular 
to  the  pipeline  axis 

(2)  Wave  period,  which  represents  the  duration  of  the  flow  in 
one  direction 

(3)  Length  of  the  horizontal  excursions  of  the  water  particles 
perpendicular  to  the  pipeline  axis  (this  quantity  is  di- 
rectly proportional  to  the  product  of  the  above  two  param- 
eters) 

Assuming  that  only  water  with  a limited  range  of  temperature  is  being 
dealt  with,  the  fluid  properties  will  be  ignored  for  the  present.  The 
orientation  angle  of  the  pipeline  can  be  handled  as  discussed  above, 
considering  only  the  components  of  the  horizontal  fluid  motions 
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perpendicular  to  the  pipeline  axis.  Since  the  length  of  the  horizontal 
water  particle  excursions  is  directly  proportional  to  the  product  of 
the  wave  period  and  the  maximum  horizontal  water  particle  velocity,  cnly 
four  independent  variables  are  left:  diameter,  clearance,  horizontal 
water  particle  velocity,  and  wave  period.  Thus,  any  single  parameter 
used  to  relate  C^,  (j>,  or  k to  the  wave  and  pipeline  conditions  must 
include  these  four  variables.  This  constraint  is  necessary  if  the  rela- 
tionship is  expected  to  be  valid  for  general  application  under  any  set 
of  wave  and  pipeline  conditions . 

The  four  variables  can  be  arranged  into  several  dimensionless  param- 
eters. The  important  parameters  should  include  the  following: 

(1)  relative  clearance,  clear/Dia 

where  clear  = bottom  clearance 
Dia  = pipe  diameter 

(2)  Keulegan-Carpenter  parameter,  ufflax  T/Dia 

where  T » wave  period 

umax  = comPonent  maximum  horizontal  water  particle 
velocity  perpendicular  to  the  pipeline  axis 

(3)  clear/ufflax  T 

NOTE. --Not  all  of  these  parameters  are  necessary  to  describe  the  system 
since  some  are  redundant,  but  some  may  be  more  useful  than  others. 

Since  viscosity  is  an  important  variable  involved  in  the  choking 
phenomenon,  the  Reynolds  number,  Dia/v,  and  a Reynolds  number  for 

the  clearance,  u^ax  clear/v,  are  also  important  parameters  (where  v = 
kinematic  viscosity) . 

The  dimensionless  parameters,  clear/u,,,^  T,  u^.^  T/Dia,  Uj,^  clear/v, 
and  Umax  Dia/v,  were  plotted  versus  the  lift  force  parameters,  C^,  <j>,  k, 
Cl(1-k)«  and  C^(k),  for  constant  values  of  the  relative  clearance, 
clear/Dia.  The  correlation  was  not  good  with  the  parameters  involving 
the  coefficient  of  lift  (C^,  C^Q-k),  and  C^Ck)).  However,  good  cor- 
relation was  found  between  several  of  the  dimensionless  parameters  and 
the  quantities  <ji  and  k. 

The  parameter,  clear/Ujjj^  T,  exhibited  the  best  correlation  with 
both  <j>  and  k for  each  relative  clearance,  although  there  was  some  varia- 
tion in  these  relationships  for  the  data  corresponding  to  the  different 
pipe  diameters  (see  Figs.  47  to  52).  Although  the  differences  are  not 
large,  tuc  data  do  indicate  the  presence  of  a scale  effect  in  these 
relationships . 
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<$>  and  k were  also  correlated  with  the  Keulegan-Carpenter  parameter, 
%ax  T/Dia.  However,  these  relationships  were  not  the  same  when  the 
data  corresponding  to  a given  relative  clearance  were  compared  for  dif- 
ferent pipe  diameters.  The  relationships  were  the  same  for  a given 
absolute  clearance,  rather  than  a relative  clearance  (clear/Dia) . These 
relationships  are  shown  in  Figures  53  and  54  for  the  combined  data  from 
all  three  pipe  diameters. 

The  parameter,  Uj^  clear/v,  demonstrated  correlation  with  both  <{» 
and  k,  but  these  relationships  also  exhibited  a scale  effect,  such  that 
the  relationships  for  a given  relative  clearance  were  not  the  same  when 
comparing  the  data  for  different  pipe  diameters . Figures  55  and  56  are 


examples  of  these  relationships  for  the  4-inch-diameter  pipeline. 

Correlation  between  the  Reynolds  number,  Dia/v,  and  the  param- 

eters, and  k,  was  not  good,  especially  when  comparing  the  data  fir 
the  different  pipe  diameters. 


Since  none  of  the  above  dimensionless  parameters  alone  could  be 
used  to  determine  a value  of  (J)  or  k for  any  given  pipe  diameter,  clear- 
ance, and  wave  condition  due  to  the  presence  of  scale  effects,  several- 
of  the  parameters  were  combined  in  various  ways  to  form  different  dimen- 
sionless parameters  containing  all  four  of  the  important  variables 
(clear,  L)ia,  u^*,  and  T) . An  attempt  was  made  to  find  a single  param- 
eter containing  all  of  the  important  variables  that  was  well  correlated 
with  <♦>  or  k for  all  wave  conditions,  pipeline  sizes,  and  configurations. 

Several  relationships  were  found  that  exhibited  good  correlation  for 
all  the  wave  and  pipeline  conditions  tested.  However,  since  this  is  a 
model  study  and,  therefore,  limited  to  lower  values  of  the  Keulegan- 
Carpenter  parameter  and  Reynolds  number  than  prototype  design  situations 
in  the  ocean,  caution  should  be  used  in  extrapolating  these  results. 

The  dimensionless  combination,  (clear/u_ax  T) (Dia/umaxT) , demon- 
strated the  best  correlation  with  both  <|>  ana  k for  all  conditions 
tested.  These  relationships  are  given  in  Figures  57  and  i>8 . Since  both 
k and  (j)  define  the  point  at  which  choking  occurs  in  the  wave  cycle,  it 
appears  that  the  choking  phenomenon  is  directly  dependent,  on  the  water 
particle  excursions  relative  to  both  the  pipe  diameter,  (Dia/u^y  T) , 
and  the  bottom  clearance,  (clear/ T) . 

Although  the  parameter,  (clear/u^y  T),  is  equivalent,  to  the  ratio 
of  the  bottom  clearance  to  the  horizontal  excursion  of  the  water  parti- 
cles (differing  only  by  the  constant  1/tt)  , the  quantity  (u^y  T)  should 
not  be  thought  of  only  as  defining  the  length  of  the  water  particle 
excursions.  Both  variables,  u,,,^  and  T,  are  independently  important  in 
defining  the  choking  phenomenon.  The  larger  Umax*  socner  the  chok- 
ing conditions  will  develop  in  the  wave  cycle  for  a given  clearance  and 
pipe  diameter.  Similarly,  since  the  wave  period,  T,  defines  the  duration 
of  the  horizontal  flow  in  one  direction,  the  larger  the  wave  period,  the 
sooner  choking  will  develop  relative  to  the  temporal  length  of  the  wave 
cycle. 
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Figure  55.  k versus  ■— a*~- for  4-inch  diameter 
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The  slight  amount  of  scatter  in  these  plots  in  the  vicinity  of  k » 

1 and  <j>  ■ 90°  is  due  to  the  error  in  calculated  values  of  <J>  and  k for 
the  largest  bottom  clearances  where  the  lift  effect  was  small  (as  dis- 
cussed above) . 

Larger  values  of ‘the  dimensionless  combination,  (clear/uroax  T) 

(Dia/ii^py  T),  than  given  in  the  plots  would  correspond  to  larger  bottom 
clearances  and  pipe  diameters  relative  to  the  maximum  velocities,  wave 
periods,  and  water  particle  excursions.  For  these  conditions,  the 
values  of  k and  $ would  remain  at  1 and  90°,  respectively,  while  the 
lift  effect  would  eventually  diminish  to  zero  with  increasing  values  of 
this  parameter.  These  trends  are  evident  in  the  data  token  at  the 
largest  bottom  clearances  (1  and  2 inches),  although  these  data  were 
not  included  in  the  ?bove  plots. 

Similarly,  lower  values  of  the  dimensionless  parameter  than  given 
in  the  plots  would  correspond  to  higher  maximum  velocities,  wave  peri- 
ods, and  water  particle  excursions  relative  to  the  smallest  bottom 
clearances  and  pipe  diameters.  So  for  lower  values  of  this  parameter, 
both  k and  $ should  remain  at  their  defined  minimum  values  of  0 and  0°, 
respectively,  corresponding  to  lift  forces  acting  in  the  upward  direc- 
tion only,  with  very  little  or  no  flow  possible  under  the  pipe  section. 

Although  $ was  defined  as  varying  from  D*  to  90°  only,  negative 
values  of  <j>  are  exhibited  in  the  data  for  the  lowest  values  of  the 
dimensionless  parameters  plotted.  However,  since  most  of  these  data 
points  correspond  to  the  smallest  diameter  pipeline  model  tested  (2 
inches),  this  could  be  partly  due  to  experimental  error,  since  the  measured 
forces  *ere  smallest  for  the  smallest  model.  Also,  part  of  this  dis- 
crepancy could  be  due  to  the  difficulty  of  accurately  defining  the  pe.'k 
of  the  wave  crest  in  the  experimental  wave  records.  This  point  was 
arbitrarily  defined  as  the  midpoint  of  the  zero  crossings  on  either 
side  of  the  wave  crest  in  the  digitized  data  records.  However,  in  some 
cases,  the  waves  were  not  perfectly  symmetrical,  so  the  maximum  eleva- 
tion of  the  water  surface  did  not  coincide  exactly  with  the  midpoint  of 
the  zero  crossings.  This  was  especially  true  of  the  largest  waves  with 
the  longest  periods,  which  in  the  plotted  relationships  would  correspond 
to  the  minimum  values  of  the  dimensionless  parameters  (at  the  lowest 
bottom  clearance  tested).  Thus,  the  actual  kinematics  under  these  waves 
would  be  slightly  out  of  phase  with  the  calculated  kinematics,  resulting 
in  an  error  in  the  calculated  value  of  <fr.  However,  this  source  of  error 
should  be  the  same  for  the  large- diameter  models  as  for  the  smallest  models. 

5 . Relationships  Between  <j>  (clear/Dia)  and  k (clear/Dia)  and  Parameters 

Defining  the  Wave  and  Pipeline  Conditions. 

Many  other  useful  relationships  were  found  by  multiplying  <p  and  k 
by  the  relative  clearance,  (clear/Dia),  and  plotting  these  dimensionless 
products  versus  various  dimensionless  parameters  defining  the  wave  and 
pipe  conditions.  Figures  59  to  62  are  examples,  although  several  other 
parameters  also  showed  good  correlation  with  (ft  (clear/Dia)  and  k (clear/ 

Dia) . 
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Both  <j>  (clear/Dia)  and  k (clear/Dia)  are  correlated  with  the  dimen- 
sionless combinations  (clear/u^,^  T)  (clear/Dia)  and  /Dia/umax  T 
(clear/Uj,  T)  (clear/Dia)  . However,  k (clear/Dia)  appears  to  be  better 
correlated  with  the  first  parameter,  while  <j)  (clear/Dia)  shows  better 
correlation  with  the  second  parameter. 

It  is  clear  that  for  values  of  the  dimensionless  parameters  lower 
than  those  shown  on  the  plots,  both  (clear/Dia)  and  k (clear/Dia) 
will  remain  at  a value  of  zero.  This  would  correspond  to  situations 
where  the  clearance  was  minimal  relative  to  the  hori.  ontal  velocities, 
wave  periods,  and  horizontal  excursions  of  the  water  particles.  Thus, 
both  k and  <j>  would  be  expected  to  equal  zero  and  0°,  r jspectively,  and 
the  relative  clearance  would  either  equal  or  approach  zero. 

Large  values  of  the  dimensionless  parameters  correspond  to  situa- 
tions where  the  clearance  is  large  relative  to  the  horizontal  veloci- 
ties, wave  periods,  and  horizontal  excursions  of  the  water  particles. 

For  these  cases,  k and  <J>  will  remain  at  maximum  values  of  1 and  90°, 
respectively,  while  the  relative  clearance,  (clear/Dia),  will  increase 
with  increasing  values  of  the  dimensionless  parameters.  But  as  the 
relative  clearance  is  increased  beyond  this  point,  the  lift  forces  will 
decrease  to  zero,  so  extension  of  the  plotted  relationships  to  much 
larger  values  of  the  dimensionless  parameters  is  of  little  value. 

6 . Relationships  Between  the  Coefficients  of  Lift  and  Parameters 

Defining  the  Wave  and  Pipeline  Conditions . 

The  coefficient  of  lift,  C^,  the  effective  positive  coefficient  of 
lift,  (^(l-k),  the  effective  negative  coefficient  of  lift,  Cj(k),  and 
the  maximum  effective  coefficient  of  lift  (maximum  of  C^fl-k)  or  Cg(k)) 
were  plotted  against  various  combinations  of  the  dimensionless  param- 
eters. The  parameter,  (clear/Un,^  T) (Dia/umaxT)  , which  previously 
gave  the  best  correlations  with  5 and  k also  demonstrated  the  best  cor- 
relation with  C^,  Ck(l-k),  and  (^(k).  However,  these  relationships 
exhibited  more  scatter  than  the  previously  discussed  interrelationships 
between  the  coefficients  of  lift  and  the  parameters,  k and  <J>,  so  it  is 
suggested  that  the  previously  discussed  relationships  be  used  for  design 
purposes . 

7.  Relationships  Between  the  Lift  Forces  and  Parameters  Defining  the 

Wave  and  Pipeline  Conditions. 

As  with  the  coefficient  of  lift,  the  total  lift  force  (FL  = 

1/2  C,  p A %ax2)  can  be  partitioned  into  the  maximum  positive  lift, 
FL(l-k),  and  the  maximum  negative  lift,  FL(k)  (Fig.  6).  These  three 
forces,  as  well  as  the  maximum  lift  force  (maximum  of  either  F^(l-k)  or 
FL(k))  were  plotted  against  various  combinations  of  the  dimensionless 
parameters.  Only  one  relationship  exhibited  good  correlations  for  the 
data  from  all  three  diameters  plotted  together.  This  was  the  Reynolds 
number,  u,j,axDia/v,  versus  the  maximum  lift  force  (either  FL(l-k)  or  FL(k) 
whichever  is  greater)  (Fig.  63). 
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Figure  63.  Maximum  lift  force  (positive  or  negative)  versus  the  Reynolds  number. 


This  relationship  shows  that  for  any  pipe  diameter,  orientation  angle, 
or  bottom  clearance,  the  maximum  lift  force  increases  with  the  Reynolds 
number  in  a regular  manner,  at  least  over  the  range  of  the  data  in  this 
investigation.  The  maximum  lift  force  may  occur  in  either  the  upward  or 
downward  direction,  depending  on  the  magnitude  of  the  bottom  clearance 
relative  to  the  wave  conditions  and  pipe  size.  This  relationship  does 
not  hold  for  the  maximum  upward  lift  or  maximum  downward  lift  alone,  but 
only  for  the  largest  of  these  two  forces  in  any  given  situation. 

8 . Relationships  Involving  the  Vertical  Coefficients  of  Mar. s anu  Drag 

and  the  Vertical  Inertial  and  Drag  Forces. 

Both  the  vertical  coefficient  of  mass  and  the  vertical  inertial 
forces  were  plotted  against  several  dimensionless  parameters  defining 
the  wave  and  pipeline  conditions,  but  no  useful  relationships  were 
found.  This  is  not  surprising  when  considering  that  the  vertical  iner- 
tial forces  are  relatively  small,  and  thus  subject  to  error  from  the 
transverse  eddy- induced  forces  which  were  not  accounted  for  in  the 
least  squares  analysis. 

No  attempt  was  made  to  plot  relationships  involving  the  vertical  drag 
forces  or  drag  coefficients,  since  these  forces  were  negligible. 

9 • Relationships  Between  the  Horizontal  Coefficient  of  Mass  and 

Parameters  Describing  the  Wave  and  Pipeline  Conditions. 

A limited  number  of  horizontal  force  data  were  taken  using  the  4-inch- 
diameter  two-dimensional  model.  Values  of  CM  and  CD  were  calculated 
from  the  least  squares  analysis,  and  an  attempt  was  made  to  relate  these 
coefficients  to  various  dimensionless  parameters  describing  the  wave  and 
pipeline  conditions. 

Figure  64  shows  the  horizontal  coefficient  of  mass  plotted  versus  the 
relative  clearance,  clear/Dia,  together  with  the  potential  flow  solution 
for  a circular  cylinder  in  the  vicinity  of  a plane  wall  subject  to  a 
uniform  flow  with  constant  acceleration  (Grace,  1974)  . The  data  follow 
the  potential  flow  solution  reasonably  well,  although  for  a given  rela- 
tive clearance,  there  appears  to  be  some  variation  in  the  value  of  CM 
with  varying  wave  conditions.  Also,  the  wave  force  data  give  slightly 
higher  values  of  the  coefficient  of  mass  for  the  highest  bottom  clearances 
tested.  Although  the  experimental  data  are  limited,  they  indicate  that 
the  potential  flow  solution  may  be  very  useful  in  determining  a value  for 
the  horizontal  coefficient  of  mass,  at  least  for  wave  conditions  where  the 
inertial  forces  predominate  over  the  drag  forces. 

However,  since  there  was  some  variation  in  the  values  of  CM  for 
different  wave  conditions  for  the  same  relative  clearance,  an  attempt 
was  made  to  determine  relationships  between  the  horizontal  coefficient 
of  mass  and  the  various  dimensionless  parameters  defining  the  wave  and 
pipeline  conditions.  Reasonably  good  correlations  were  found  between 


several  of  the  parameters.  Figure  65  shows  the  relationship  for  Cm 
versus  clear/umaxT. 

10.  Relationships  Involving  the  Horizontal  Coefficient  of  Drag. 

The  horizontal  coefficient  of  drag  was  plotted  against  several 
dimensionless  parameters,  but  no  useful  relationships  were  found.  This 
was  expected  since  the  horizontal  drag  forces  in  this  investigation  w,jre 
much  smaller  than  the  inertial  torces,  due  to  the  horizontal 

excursions  of  the  water  particles  relative  to  the  diameter  of  the  pipe- 
line. 


11.  Example  Problems. 


GIVEN:  A design  wave  with  height,  H = 10  feet  and  period,  T = 10  seconds 
acts  on  a pipeline  with  a diameter,  Dia  = 8 feet  in  a water 
depth,  d = 80  feet.  The  pipeline  is  oriented  at  an  angle  of 
30°  with  respect  to  the  wave  crests.  Section  A of  thf  peline 
is  in  contact  with  the  bottom;  section  B spans  the  bottom  at  a 
clearance,  clear  = 6 inches. 


FIND:  For  both  sections  A and  B,  find 

(a)  the  values  of  the  lift  force  parameters  (CL,  <fr,  and  k); 

(b)  the  maximum  positive  and  negative  lift  forces; 

(c)  the  positions  of  these  maximum  lift  forces  in  the  wave 
cycle;  and 

(tl)  the  lift  force  at  8 = 120°  in  the  wave  cycle. 


SOLUTION : 

L, 


m*Jlu 

2 TT 


5.12  (10)2  = 512  feet 


80 

512 


0.1562 


j on 

Using  tables,  j-  = 0.1885,  so  L = - ■ ^ = 424  feet 
sinh  ~ = 1.481 

z = distance  from  bottom  to  center  of  pipe 
sections. 


i 
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Figure  65.  Horizontal 


yprtSHvr:  »* 


For  section  A (clear  = 0) 
z = 4 feet 


L " 424 


= 0.00943 


From  tables,  cosh 


1.0017 


r ! - 

i v( 

i\  \ 


ir  H COshC  LZ)  tt(10) (1.0017)  , ,,  _ _ 

umax  ~ T s inh (— — — )~  ' (WTO!!)  * 2,12  f t per  Second 

L 

Component  of  perpendicular  to  the  pipeline  axis  is 
Umay  (cos  30°)  « (2 . 12) (0 . 866)  = 1.84  feet  per  second 

(a)  Since  the  pipe  is  in  contact  with  the  bottom,  (clear  « 0) , 

<J>  ® 0°  and  k = 0.  From  Figure  40,  Cl  * 4.5. 

(b)  Maximum  positive  lift  (per  ur.it  length) 

FL(l-k)  - | CL  p A umax2  (1-k) 

« | (4.5)  (2) (8) (1.84)2  (1-0) 

* 121.9  pounds  per  foot. 

Maximum  negative  lift  (per  unit  length) 

Since  k = 0,  there  is  no  negative  lift,  and  the  lift  force  is 
positive  throughout  the  wave  cycle. 

(c)  Since  <j>  = 0°,  the  positive  lift  forces  are  maximum  at  0°  and  180° 
in  the  wave  cycle  (under  the  crests  and  troughs) , corresponding  to 
the  points  of  maximum  horizontal  velocities. 

The  lift  does  not  become  negative,  but  diminishes  to  zero  at  90° 
and  2708,  the  positions  of  horizontal  flow  reversal  in  the  wave 
cycle. 

(d)  At  6 = 120° 

FL  = j CL  p A u^  [cos2  (6  - <j>)  - k] 

= ~ (4.5)  (2)  (8)  (1 . 84)2  [cos2  (120°  - 0°)  - 0] 
s 30.5  pounds  per  foot 


For  section  B (clear  t;  6 inches) 
z * 4.5  feet 


z_ 

L 


4.5 

424 


0.0106 


From  tables,  cosh 


2tti 


1.0022 


. ,2irz. 

, H “shf-r) 

“«”**  ‘ T ci„v,2ltd. 

sinh(— j—) 


rTs'n1  85  2.13  ^eet  per  secon<^ 


(10)  (1.481) 


► ’ 

clear 

Dia 

u T 

u T 

max 

_ max 

component  of  ujnax  perpendicular  to  the  pipeline  axis  is 
%ax  Ceos  30°)  * (2. 13)  (0.866)  » 1.84  feet  per  second 
(a)  Use  Figure  57  to  determine  a value  for  k 

C9-‘  51.LILL _ o 0118 

(1.84^(10) (1.84) (10)  u.uiiB 

so  from  Figure  57,  k * 0.67 
and  from  Figure  58,  $ - 45°. 

Alternatively,  either  $ or  k could  be  determined  from  Fig.  39, 
once  the  other  is  known. 

From  Figure  46,  for  k = 0.67, 


so  CL 


CL(l-k)  = 2.75 
2.75 


(1  - 0.67) 


8.3. 


(b)  Maximum  positive  lift  (per  unit  length) 
FL(l-k)  - $ CL  p A Vx!  (1-k) 


= ~ (8. 3) (2) (8) (1. 84)2  (1  - 0.67) 


74.2  pounds  per  foot 


Maximum  negative  lift  (per  unit  length) 
" ^(k)  * “ j Cl  P A umax2  (k) 

* - j (8. 3) (2) (8) (1.84) 2 (0.67) 


150.6  pount  per  foot 


(c)  Since  0 ■ 45° , the  positive  lift  forces  are  maximum  at  0*  ♦ 45°  » 
45*  and  180°  * 45°  * 225°  in  the  wave  cycle,  and  the  negative  lift 
forces  are  maximum  at  90°  + 45°  = 135°  and  270°  + 45°  * 315°  in 
the  wave  cycle. 


(d) 


At  8 * 120° 

FL  5 | CL  p A u^,2  [cos2  (120s  - 45°)  - 0.67] 


- J (8. 3) (2) (8) (1.84)2 


[cos2  (120°  - 45°)  - 0.67] 


= - 135.6  pounds  per  foot 

Again,  it  should  be  stressed  that  the  relationships  involving  the 
lift  force  parameters,  Cl,  <J>,  and  k,  were  determined  from  model  studies 
conducted  r.t  much  lower  values  of  the  Keulcgan- Carpenter  parameter  and 
Reynolds  number  than  those  encountered  in  full-scale  situations  in  the 
ocean  Therefore,  caution  should  be  used  in  extrapolating  these  results 
to  prototype  designs. 

Further  studies  using  a larger  scale  facility  are  necessary  to 
evaluate  the  importance  of  rcale  effects  in  these  relationships,  to 
determine  their  limitations,  ai  d possibly  to  extend  or  modify  them  so 
they  are  valid  for  any  scale. 

IV,  CONCLUSIONS 

1.  The  traditional  steady- flow  lift  f rce  model,  expressed  as 
Fl  = 1/2  Cl  P A i2,  is  not  a suitable  model  for  the  description  of  wave- 
induced  lift  forces.  This  model  assumes  that  the  lift  force  acts  in 
one  direction  only  (upward  or  downward)  throughout  the  entire  wave 
cycle. 


2.  For  pif  “lines  located  at  a small  clearance  above  the  bottom, 
a viscous  choking  effect  limits  the  maximum  velocities  through  the 
constriction  formed  by  the  bottom  clearance.  Correspondingly,  the 
pressure  drop  on  the  bottom  side  of  the  pipe  section  is  also  limited. 

In  contrast,  the  flow  velocities  and  corresponding  pressure  drop 
over  the  top  side  cf  the  pipeline  are  not  limited.  As  the  choking 
effect  develops  and  the  flow  becomes  restricted  through  the  bottom 
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clearance  constriction,  more  of  the  flow  must  be  diverted  over  the  top 
of  the  pipe  section,  resulting  in  a downward  shift  in  the  stagnation 
point,  as  well  as  an  increase  in  the  flow  velocities  and  associated 
pressure  drop  cv,r  tl.s,  top  side  of  the  pipeline. 

The  induced  changes  in  the  flow  pattern,  velocities,  and  associated 
pressure  distribution  over  the  pipe  section  due  to  choking  through  the 
bottom  clearance  constriction  result  in  an  upward  lift  force,  rather 
than  the  downward  lift  force  predicted  by  potential  flow  theory. 

3.  Thus,  for  an  oscillatory  wave-induced  flow,  the  lift  force  acts 
downward  in  t.hoso  parts  of  the  wave  cycle  where  the  horizontal  water 
pa  'tide  velocities  are  not  high  enough  to  produce  choking  through  the 
bottom  clearance.  In  this  case,  the  unrestricted  flow  is  faster  through 
the  bottom  clearance  constriction  than  over  the  top  of  the  pipe  section, 
so  the  corresponding  pressure  distribution  results  in  a negative  lift 
toward  the  bottom  boundary. 

However,  in  those  parts  of  the  wave  cycle  where  the  horizontal 
velocities  are  sufficient  to  induce  choking  through  the  bottom  clearance 
constriction,  the  lift  force  acts  in  an  upward  direction. 

4.  For  a given  pipe  diameter  and  wave  condition,  as  the  bottom 
clearance  is  increased,  higher  velocities  are  necessary  to  produce  the 
choking  effect.  Thus,  the  negative  lift  force  can  reach  a greater 
magnitude  and  occur  later  into  the  wave  cycle  before  the  choking  condi- 
tion is  induced. 

Correspondingly,  the  positive  lift  that  occurs  only  after  the 
choking  condition  develops  is  limited  to  a smaller  part  of  the  wave 
cycle,  and  the  maximum  magnitude  of  these  forces  decreases  with 
increasing  clearance.  In  addition,  since  there  is  a small  timelag 
involved  in  the  development  of  the  choking  phenomenon  and  the  transition 
from  negative  to  positive  lift,  the  maximum  positive  lift  occurs  later 
into  the  wave  cycle,  although  its  magnitude  is  diminishing. 

5.  All  major  features  of  the  wave-induced  lift  force  phenomenon 
can  be  described  adequately  by  a modified  lift  force  equation,  Fl  = 

1/2  Cl  P A [cos2  (0  - $)  - k],  where  $ represents  a phase  shift 

in  the  position  of  the  maximum  positive  (upward)  lift  force  i-elative 
to  the  point  of  maximum  horizontal  velocity  at  the  center  of  the  wave 
crest,  and  k represents  the  proportion  of  the  total  lift  force  cycle 
that  acts  in  the  negative  (downward)  direction.  The  values  of  $ and 

k vary  from  0*  and  0,  respectively,  for  the  case  of  a pipeline  touching 
the  bottom,  and  increase  with  increasing  clearance  (for  a given  pipe- 
line and  wave  condition)  to  maximum  values  of  90°  and  1,  respectively, 
when  the  pipeline  is  far  enough  from  the  bottom  so  that  the  choking  condi 
tion  does  not  develop.  $ = 0°  and  k » 0 correspond  to  lift  forces  that 
are  positive  throughout  the  wave  cycle,  with  maximums  occurring  at  the 
points  of  maximum  horizontal  velocity  under  the  wave  crests  and  troughs. 

= 90*  and  k * 1 correspond  to  negative  lift  forces  throughout  the  wave 


\ cycle,  ;ith  maximum  downward  forces  occurring  under  the  crests  and 

( troughs  of  the  passing  waves.  These  two  cases  represent  the  extreme 

conditions  bounding  the  lift  force  phenomena.  At  any  intermediate 
clearance  between  these  limiting  cases,  both  positive  and  negative 
*ift  forces  will  occur  at  different  parts  of  the  wave  cycle,  and  the 
positions  of  the  maximum  upward  and  downward  lift  forces  will  not 
coincide  with  the  positions  of  maximum  horizontal  velocities  in  the 
wave  cycle. 

In  order  to  use  this  lift  force  model,  values  of  the  parameters, 
i . C^,  and  k,  must  be  determined  for  the  given  set  of  wave  and  plpe- 

j line  conditions.  A model  investigation  was  carried  out  to  determine 

relationships  between  these  parameters  and  various  dimensionless  param- 
( # eters  defining  the  wave  and  pipeline  conditions. 

| 6.  A direct  relationship  was  found  between  the  lift  force 

parameters,  $ and  k.  Relationships  were  also  found  between  the 
I coefficient  of  lift,  C^,  and  both  $ and  k.  In  addition,  CL  can  be 

• partitioned  into  the  positive  effective  coefficient  of  lift,  C.  (1-k), 

{ and  the  effective  negative  coefficient  of  lift,  C^fk).  Both  of  these 

parameters  are  also  related  to  both  <£  and  k.  The  correlation  is  better 
with  k than  <J>  for  the  relationships  involving  CL,  CL(l-k),  and  C^k). 

All  of  these  relationships  were  the  same  for  all  pipe  diameters, 

! bottom  clearances,  and  wave  conditions  tested. 

i 7.  The  average  value  of  CL  at  k ■ 0 and  <f>  » 0“  (which  corresponds 

j to  a pipeline  in  contact  with  the  bottom  with  no  clearance)  is  4.5. 

| This  is  the  same  as  the  potential  flow  solution  for  the  lift  force  on 

\ a circular  cylinder  against  a plane  wall  subject  to  a steady,  inviscid 

I flow  parallel  to  the  wall. 

8.  Maximum  values  of  CL  occur  at  k = 1/2  and  d>  = 30°,  where 
Cl  * 9.  In  the  interval  from  k * 0 to  1/2  and  <J>  = 0°  to  30°,  the 
j effective  positive  coefficient  of  lift  CL(l-k)  remains  at  approximately 

■ 4.5,  while  the  effective  negative  coefficient  of  lift  C, (k)  increases 

from  0 to  4.5.  In  the  interval  from  k = 1/2  to  1 and  <f>  = 30°  to  90°, 
CL(l-k)  decreases  to  0,  while  Cj^ Ck)  increases  to  reach  a maximum  of 
about  6 or  7 at  k * 0.75  and  <j>  = 45°,  and  then  decreases  to  a maximum  of 
4.5  at  k = 1 and  (j)  - 90#. 

\ 

, ‘ 9.  Using  the  above  relationships  between  Cl,  <j>,  and  k,  if  either 

♦ or  k is  known,  the  remaining  two  parameters  can  be  determined. 
Therefore,  an  attenpt  was  made  to  find  relationships  between  $ and  k 
and  various  dimensionless  parameters  defining  the  wave  and  pipeline # 
conditions. 

The  best  correlation  was  found  in  the  relationships  between  $ and 
k and  the  parameter  clear/ u^^T  for  constant  values  of  the  relative 
clearance,  clear/Dia.  However,  comparison  of  the  data  corresponding  to 
the  different  pipe  diameters  indicates  a slight  scale  effect  is  present. 
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♦ and  k were  also  related  to  the  parameter  u^Mclear/v  for  constant 
values  of  clear/Di a,  although  the  scale  effect  was  worse  for  these  rela- 
tionships. <j>  and  k showed  very  good  correlation  with  the  Keulegan- 
Carpenter  parameter,  ul  T/Dia,  although  these  relationships  were  the 
same  for  a constant  absolute  clearance,  rather  than  a constant  relative 
clearance.  Correlation  between  $ and  k and  the  Reynolds  number  was 
poor. 


10.  Because  a scale  effect  was  evident  in  the  above  relationships, 
several  of  the  dimensionless  parameters  were  combined  to  form  new 
dimensionless  parameters  that  contained  all  of  the  ing>ortant  variables 
(clear,  Dia,  u^x*  and  T) . An  attempt  was  made  to  find  a single  param 
eter  that  was  related  to  either  ^ or  k for  all  wave  and  pipeline  condi- 
tions tested  in  this  investigation. 


Both  $ and  k showed  very  good  correlation  with  the  parameter 
(clear/u^T)  (Dia/u-^T) . These  relationships  were  valid  for  all  pipe 
diameters;  bottom  cTearam 
tested. 


Learances,  orientation  angles,  and  wave  conditions 


In  addition,  the  relative  clearance  was  combined  with  both  and 
k to  form  the  quantities  $(clear/Dia)  and  k (clear/Dia) , both  of  which 
exhibited  very  good  correlation  with  more  of  the  dimensionless  combina- 
tions than  either  $ or  k alone.  k(clear/Dia)  was  best  correlated  with 
(clear/u^xT)  (clear/Dia) . 4>(clear/Dia)  was  also  correlated  with  this 
parameter,  but  exhibited  better  correlation  with  the  parameter 
✓bia/umayT  (clear/u^T)  (clear/Dia) . 

11.  CL,  CL(l-k),  ana  C^(k)  were  correlated  with  the  same  parameter 
as  $ and  k,  (clear/i^^T)  (Dia/v^  T) , However,  these  correlations  were 
not  as  good  as  the  previous  correlations  between  the  coefficients  of 
lift  and  k or  <j>. 

12.  For  a pipeline  that  is  not  parallel  to  the  wave  crests,  the 
lift  forces  are  apparently  due  only  to  the  components  of  the  horizontal 
water  particle  velocities  perpendicular  to  the  axis  of  the  pipeline. 
Using  this  convention,  consistent  values  of  the  coefficient  of  lift, 

Cj,,  are  obtained  for  all  angles  of  orientation.  In  addition,  the  rela- 
tionships between  the  lift  force  parameters  Cl,  $,  and  k,  as  well  as 
relationships  between  these  parameters  and  various  dimensionless  param- 
eters defining  the  wave  and  pipeline  conditions,  are  identical  for  all 
angles  of  orientation. 

13.  The  maximum  lift  force  (FL(l-k)  or  FL(k) , whichever  is  greater) 
exhibited  good  correlation  with  the  Reynolds  number,  (umaxD*a/v) • This 
relationship  did  not  hold  for  the  maximum  positive  lift  tFr (l-k) ) or 
the  maximum  negative  lift  (FL(k))  alone,  but  only  for  the  largest  of 
these  two  forces  in  any  situation.  The  relationship  was  the  same  for 
all  diameters  over  the  range  of  conditions  tested. 


(20 


a 


14.  The  horizontal  coefficient  of  mass,  C^t  showed  excellent 
agreement  with  the  potential  flow  solution  for  a circular  cylinder  in 
the  vicinity  of  a plane  wall  subject  to  a uniform  flow  with  constant 
acceleration.  These  results  indicate  that  the  potential  flow  solution 
may  be  useful  for  selecting  a value  of  for  wave-induced  forces,  at 
least  for  situations  in  which  the  inertial  forces  predominate  over  the 
drag  forces.  The  horizontal  C|^  was  also  correlated  with  several  of  the 
dimensionless  parameters  defining  the  wave  and  pipeline  conditions, 
such  as  the  parameter  clear/t^^T. 

V.  RECOMMENDATIONS  FOR  FURTHER  RESEARCH 

1.  Experiments  similar  to  this  investigation  should  be  carried  out 
in  a larger  wave  tank  facility.  This  would  allow  the  testing  of  larger 
diameter  pipeline  models  as  well  as  experiments  at  higher  Reynolds 
numbers  and  higher  values  of  the  Keulegan- Carpenter  parameter.  Such 

an  investigation  is  necessary  to  determine  the  validity  of  extrapolat- 
ing the  results  of  the  present  study  to  design  situations  in  the  ocean, 
and  to  point  out  any  weaknesses  or  limitations  of  the  proposed  lift 
force  model  due  to  scale  effects. 

2.  It  would  be  of  interest  to  perform  experiments  to  evaluate  the 
magnitude,  phase,  and  frequency  spectra  of  the  vertical  transverse 
lift  forces  due  to  eddy  shedding  for  a horizontal  cylinder  subject  to 
oscillatory  horizontal  flow  velocities.  This  could  be  done  by  oscil- 
lating a test  cylinder  horizontally  in  still  water  away  from  a boundary, 
or  by  using  a pulsating  flume  facility.  The  horizontal  flow  patterns 

at  the  bottom  could  be  simulated,  but  without  the  lift  force  phenomenon 
due  to  the  boundary.  Only  the  transverse  lift  forces  due  to  eddy 
shedding  would  act  in  the  vertical  direction,  so  the  magnitude  and  time 
history  of  these  forces  could  be  easily  measured. 

A thorough  analysis  of  the  eddy  forces  for  different  pipe  diameters 
and  flow  conditions  would  allow  an  evaluation  of  their  importance  rela- 
tive to  the  Bernoulli-type  lift  forces,  and  at  the  same  time  explain 
some  of  the  variations  in  the  vertical  wave  force  parameters  calculated 
from  an  analysis  which  neglected  the  eddy  forces  because  they  could 
not  be  separated  analytically  because  of  their  random  nature.  Adequate 
knowledge  of  the  eddy  forces  would  allow  the  addition  of  the  eddy  lift 
force  term,  * 1/2  Cj^  p A uTnax2 , to  the  Morison  equation  with  appro- 
priate values  of  the  coefficient  C[  for  any  given  set  of  wave  and  pipe- 
line conditions. 

It  should  be  noted  t.iat  evaluation  of  the  eddy  forces  for  a cylinder 
away  from  a boundary  would  only  give  an  approximate  estimate  of  the 
eddy  release  phenomenon  for  a pipe  located  near  the  bottom.  The  presence 
of  the  bottom  boundary  changes  the  flow  pattern,  velocities,  and  pressure 
distribution  around  the  cylinder,  and  therefore  would  be  expected  to 
have  some  effect  on  the  formation  and  release  of  eddies. 
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3.  Since  the  restricted  flow  through  the  narrow  bottom  clearance 
constriction  is  the  critical  part  of  tho  lift  force  phenomenon,  the 
effect  of  pipeline  roughness  and  bottom  roughness  on  the  wave -induced 
lift  forces  should  be  studied.  This  has  practical  significance,  since 
the  ocean  floor  is  not  necessarily  smooth,  and  pipelines  installed  in 
marine  waters  may  soon  becoire  encrusted  with  marine  organisms,  thus 
increasing  their  surface  roughness . 

4.  The  effect  cn  the  lift  force  phenomenon  of  a horizontal  bottom 
current  superimposed  or.  the  oscillatory  motions  of  the  wave  action 
should  be  investigated. 

5.  The  effect  of  porosity  of  the  bottom  on  lift  forces  should  also 
be  investigated. 
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LEAST  SQUARES  ANALYSIS  OF  EXPERIMENTAL  DATA 


Using  Moris  on’s  method  for  the  calculation  of  wave  forces  on  a 
pipeline,  the  vertical  component  of  the  wave-induced  force  may  be 
expressed  as: 


Fv  = (Fl)v  + (fd)v  + fL  + fL 

= cM  p V ♦ 1/2  cD  p A v !v! 

+ 1/2  CL  p A Umax2  [cos2  (6  - <p)  - k] 

+ 1/2  C[  p A umax2.  (Al) 

Since  the  transverse  lift  force  associated  with  eddy  shedding  (F^)  is 
a random  phenomenon,  there  is  no  way  to  handle  its  time  history  in 
analyzing  a wave  force  record  with  several  other  forces  occurring 
simultaneously.  Because  the  Bernoulli -type  lift  forces  were  much 
larger  than  the  eddy-associated  forces  for  a pipeline  located  close  to 
the  bottom,  the  eddy-associated  lift  force  term  was  dropped  from  the 
analysis . 

The  vertical  components  of  the  water  particle  velocities  and 
accelerations  near  the  bottom  are  small  in  comparison  with  the  corre- 
sponding horizontal  components.  As  a result,  the  vertical  lift  forces 
due  to  the  horizontal  components  of  the  water  particle  velocities  are 
generally  much  larger  than  the  vertical  drag  and  inertial  forces.  The 
drag  forces  are  especially  insignificant  since  the  vertical  excursions 
of  the  water  particles  near  the  bottom  are  smaller  than  the  diameter 
of  the  pipeline. 


(A 9) 


pATV2H2 

j_ 

2T 


,2  /27T2» 
cosh  (~l~) 

• t_  2 /-2ird, 
smh  (-j~ ) 


The  expressions,  F^v,  Fpy,  and  Fpv,  are  constant  for  a given  set  of 
wave  and  pipeline  conditions. 

Linear  wave  theory  was  used  in  the  analysis  because,  as  discussed 
previously,  there  seems  to  be  no  obvious  way  of  accurately  describing 
the  lift  force  phenomenon  mathematically  using  higher  order  theories. 

Since  the  lift  forces  much  larger  than  the  vertical  drag  or  inertial 
forces,  with  the  drag  forces  being  almost  completely  insignificant, 
there  was  no  point  in  using  higher  theories  to  eApress  the  vertical 
components  of  the  drag  and  inertial  forces. 

For  any  vertical  wave  force  record  in  which  the  corresoonding  wave 
and  pipeline  conditions  are  known,  a least  squares  analysis  can  be 
performed  on  the  data  to  determine  the  values  of  the  unknown  parameters 
Cl,,  *P>  k,  C^,  and  Cp  in  the  vertical  wave  force  equation.  The  least 
squares  analysis  yields  the  values  of  these  five  parameters  which  best 
fit  the  force  data  throughout  the  entire  wave  cycle.  This  is  accomplished 
by  determining  the  values  of  these  parameters  which  minimizes  the  sum  of 
squares  of  the  difference  between  the  observed  force  data  and  the  corre- 
sponding forces  calculated  with  the  mathematical  model  throughout  a 
complete  wave  cycle. 

Using  the  appropriate  trigonometric  identities. 


cos2  (8  - <fr)  = 1/2  + 1/2  cos  2 (0  - <j>) 

= 1/2  + 1/2  (cos  20  cos  2<j>  + sin  20  sin  2<{>) , (A10) 


so  the  lift  force  equation  can  be  expressed  as: 

Fp  = l/2CppAumax2 [l/2cos  2<j>  cos  20  + l/2sin  2<f>  sin  20  + 1/2  - k]  (All) 

or  Fp  = Aj  cos  20  + Bp  sin  20  + Cp  (A12) 

where  Ai  = 1/4  Cp  p A cos  2 <J>  = 1/2  Cp  Fpv  cos  2<p  (A13) 

Bp  = 1/4  Cp  p A Uj^2  sin  2(J>  = 1/2  Cp  Fpv  sin  2<f>  (A14) 
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2 


(1/2  - k)  = CL  FLv  (1/2  - k). 


(A15) 


C1  ' 1/2  CL  P A V 


In  an  analogous  manner,  the  vertical  components  of  the  inertial  and  drag 
forces  can  be  written  as: 


,3v 


(ft)  = c.  p (■£ ~ ) cos  e = d,  cos  e 

v rv  M ^ v3t'2max  1 


(A16) 


and  (FD)^  ■ 1/2CD  pA  vmay  sin  0 Jv^  sin  0|  = Ej  sin  0 |sin  B|,  (A17) 

3v. 


where  Dx  - CM  p (^  - - F^ 


(A18) 


(~1 

t3tJ 


. . ,2ttz. 

2itjH  sinh  t-r0 


max 


. . ,2nd. 

s inh  ("T~~) 


(A19) 


E.  = 1/2  Cn  p A v v 
1 D max  ' max1 


CD  FDv 


(A20) 


max 


. , ,2nz. 
2 m sinh 

T sinh  (~.) 


(A21) 


The  total  vertical  wave  force  at  any  position  0.  in  the  wave  cycle  can 
then  be  written  as : * 


F (0.)  = F.  + (F.)  ♦ (Fn)  = A.  cos  20.  + B,  sin  20.  + C, 

v v L v I v D'v  1 l 1 x 1 


+ Dj  cos  0i  + Ej  sin  0^^  jsin  6.J  . (A22) 


The  parameters  Aj,  Bj,  C^,  Dj,  and  Ej  are  constant  for  any  given  values 
of  C^,  <J>>  k,  C^,  and  Cq,  corresponding  to  the  particular  wave  and  pipe- 
line conditions  under  consideration. 

The  sum  of  squares  of  the  differences  between  the  observed  vertical 
forces,  F (0j.),  and  the  corresponding  calculated  forces,  Fy(0i),  is 
written  as: 
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I,  [Fv  C9t)  - Fov(0i)]’  - I [Aj  co-  26  * Bj  sin  26  ♦ C . D 
ial  i=l 

♦ EjSin  et  jsine.l  - Fov(6.)]2. 


cos  6. 


(A23) 


To  minimize  the  sum  of  squares  of  the  differences,  the  derivative  of 
this  expression  is  taken  separately  with  respect  to  each  of  the  five 
unknown  parameters  A.,  B , C , D , and  E. , and  the  resulting  expressions 
are  set  equal  to  zero,  yielding  a system1 of  five  simultaneous  equations 
with  five  unknowns.  The  system  of  equations  is  then  summed  for  each 
interval,  i,  over  a conplete  wave  cycle,  and  the  resulting  expressions 
are  solved  for  the  values  of  the  unknown  parameters  A.,  B.,  Cj,  D.,  and 
Ej  which  thus  minimize  the  sum  of  squares  or  the  differences.  The 
derivatives  are: 


SIFvfBj)  - Fqv(61)]! 


* 2A-C0S2  26.  + 2B,sin  26.  cos  26. 

1 ill  l 

♦ 2C.cos  28.  ♦ 2D, cos  8.  cos  28. 

1 i 1 i i 

+ 2E,sin8,  I sin  8.1  cos  28. 

1 i ' l 1 i 

- 2 F (6.)  cos  28.  * 0 

ovv  XJ  1 


(A24) 


a [F  (6.)  - F (8. )]2 

-v-  1 = 2A  cos  28.  sin  28.  + 2B.sin2  28. 

1 l 1 1 l 

+ 2C,sin  28.  + 2D, cos  8.  sin  28. 

1 l 1 l l 

+ 2E.sin  8.  Isin  8.1  sin  28. 

1 l 1 l ' i 

- 2F  v (8.)  sin  28.  * 0 

ov v 1 ' 1 


(A25) 


alFvtei)  - F„v(0i']; 


= 2A,cos  28.  + 2B,sin  28. 

1 i 1 i 

+ 2Cj  + 2DjCOS  8^ 

+ 2E.sin  8.  Isin  8. I 

1 l 1 i' 

- 2F  (8, ) » 0 

0‘ 


(A26) 
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W9-!1  - PovWi)ls 


* 2A.ros  20.  cos  9. 

i i i 

+ 2B1sin  20i  cos  0^^  ♦ 2C1cos  6i 

+ 2 D.cos2  0.  + 2E.sin  0.  Isin  0. | cos  0. 

1 1 1 1 1 1 

- 2F  (0.)  cos  0.  = 0 ( 

ov  1 i 


* 2A,cos  20.  sin  0.  Isin  0.1 

1 x i 1 i 

+ 2B,sin  20.  sin  0.  Isin  0.1 

1 i i ' i ' 

♦ 2C1sin  0i  | sin  Q^J 

+ 2D,cos  0.  sin  0.  Isin  0. I 

1 l l 1 l 1 

+ 2Ej  (sin  0^  | sin  0^|)2 

- 2F  (0.)  sin  0.  Isin  0-1  = 0 

ov  v l'  1 ' l1 


(A28) 


For  an  even  number  of  equally  spaced  time  intervals,  @i>  summed  over 
j.  conplete  wave  cycle,  many  of  the  terms  cancel  out  due  to  the  symmetry 
of  these  sinusoidal  functions.  Thus, 


£ cos  0.  = 0 
i=l  1 


(A29) 


y cos  20.  = 0 


(A30) 


n 

J.  sin  20.  = 0 
x=l  i 


(A31) 


r sin  0.  | sin  0. | = 0 

l i i 


(A32) 


cos  0^  cos  20^  = 0 


(A33) 
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cos  0^  sin  20^  * 0 


(A34) 


n 


l 

i*l 


n 

£ cos  0.  sin  0.  jsin  0.  | = 0 (A3S) 

i=l  1 1 


n 

J cos  20.  sin  20.  = 0 (A36) 

i=l  x 1 


n 

l cos  20.  sin  0.  jsin  0. | = 0 (A37) 

i=l  11  1 


n 

l sin  20.  sin  0.  |sin  0. | = 0 (A30) 

i=l  11  1 


when  taken  over  a complete  ways  cycle.  As  a result,  only  the  squared 
terms,  and  the  terms  involving  the  observed  forces,  F0  (0j),  remain  in 
these  equations.  The  resulting  expressions  are: 


n 

A l cos2  20. 
Ai=l  1 


n 


i=l 


F (0.) 
ov  1 


cos  20^  - 0 


(A39) 


n 

B.  7 sin2  20. 

hit 


n 

y F (0.)  sin  20.  = 0 
L,  ov  1 1 


n 


n n 

I F (0.)  = nC.-  y F (0.)  * 0 

.L.  OVv  XJ  1 .**,  OVv  l' 

1=1  1=1 


(A40) 


(A41) 


n 

n 

y cos2  0. 

- \ 

F (0.)  cos  0.  = 0 

i=l 

i=l 

OV  1 1 

(A42) 
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(A43) 


Ea  l (sin  Qi  | sin  Oj)3 


n 

I 

‘i-l 

n 


- X Fov(0i5  sin  6i  lsin  9ii  * 0 


where  n is  the  total  number  of  values  taken  from  the  vertical  wave  force 
record  (from  cn  even  number  of  equally  spaced  intervals  per  wave  cycle, 
and  over  any  number  of  complete  wave  cycles^  and  i is  the  number  of  the 
interval . 

These  expressions  are  easily  solved  for  the  unknown  paramete7  '•  A^, 
Bj,  Cy  Dj,  and  Ej,  yielding: 


n 


l CM  2ei 

i=l 


n 


J cos2  26^ 
i*l 


(A44) 
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.1  Fov">i>  3ln  29i 
1*1 
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Y sin2  29. 
i=l 


(A45) 


n 


I F (9.) 
. ovv  1 
i*l 


(A46) 


n 


y F (e.)  cos  9. 
.L.  ovv  i'  i 

i*l 


n 


E cos*  e. 
i*l 


(A47) 
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I F (8.)  sin  6 | sin  9.| 

i«l  _ 1 

l (sin  9.  jsin  9.))* 
i«l  1 


(A48) 


With  these  relationships,  the  corresponding  values  of  the  parameters 
Cj,»  <j>»  h,  C^,  and  Cjj  in  the  vertical  wave  force  equation  which  best 
fit  the  data  throughout  the  complete  wave  cycle  can  be  obtained. 

The  coefficients  of  mass  and  drag,  and  Cq,  are  obtained  directly 
from  the  parameters  rJj  and  Ej,  since 


l Fov»l>  Si 
1«1 

n 

F'  l cos2  8. 


(A49) 


CD“  -K 


! f0Vcv sin  ei  isin  eJ 

i«l 

Fdv  Jl  ^ 9i  'Sin  0i^)2 


. (A50) 


Since  Aj  * 1/2  CL  FLv  cos  2 $ and  Bj  « 1/2  FLv  sin  2$,  the  phase 
shift  parameter  $ can  be  obtained  from: 


♦ * (2*)  ’ 2/2  tan-1  - 1/2  tan'1  (ji) 

3ince  1/2  CL  F^v  cancels  out  of  the  expression  |^-j . Thus, 


(A51) 


After  <j>  is  known,  the  coefficient  of  lift,  C^,  can  be  obtained  from 
either  Aj  or  Bj,  since 


Al  l t0i)  «s  29i 

* i»l 

i*  1/2  Flv  cos  2<j)  n 

1/2  FLv  cos  2$  l cos2  20J 
i*l 


Bj  . ^ Fov  (®i^  S2n  ^®i 

CL  = ITTT^piin-^  = 1 n 

1/2  FLv  sin  24>  £ sin2  20^ 
i=i 


(AS3) 


(AS4) 


Alternatively,  could  be  obtained  from  Aj  and  Bj  directly  without  first 
solving  for  <J>,  since 


(1/2  C»  FLv  cos  2<j>)2  + (1/2  CL  FLv  sin  2<J>) 


(1/2  CL  FLv)  (cos'  24>  ♦ sin'  2(j)) 
1/2  C,  F»„  . 


n 

2 

n 

\ Fov  cos 

l Fov  (0t)  sin  20* 

i»l 

i«l 

n 

l cos2  20^ 

U-i 

♦ 

I sin2  20. 

U-l 

Finally,  the  parameter  k can  be  obtained  from  knowing  the  value 
of  Cl* .since 


k ■ 1/2  - 


CL  fLv 


1/2  - 


I Fov  CBi)  / 
i*l  / n 


CL  FLv 


(AS7) 


Thus,  once  the  vertical  wave  forces  on  a pipeline  are  measured 
experimentally,  the  values  of  the  parameters  C. , <(>,  k,  C^,  and  CD  of 
the  vertical  wave  force  equation  which  best  fit  the  data  throughout  the 
entire  wave  cycle  can  be  determined  for  the  particular  set  of  wave  and 
pipeline  conditions  tested. 

In  an  analogous  manner,  the  least  squares  analysis  car  be  applied 
to  the  horizontal  wave  force  data.  Oroitting  the  horizontal  force 
associated  with  eddy  shedding,  t’.^e  horizontal  component  of  the  wave- 
induced  force  can  be  expressed  as  equation  (2) : 


Fh  - (Fl)h  + (pD)h  * c„  p V f-  ♦ 1/2  CD  p A u|u| 


The  data  from  the  horizontal  force  measurements  show  that  the  horizontal 
eddy  forces  are  insignificant  in  comparison  to  the  horizontal  drag  and 
inertial  forces  for  the  e.^erimental  conditions  tested. 

Using  linear  wave  theory,  vhe  horizontal  components  of  the  wave 
kinematics  with  respect  to  time  can  be  expressed  as: 


tth  tjjr? 

T sinh 


cos  6 


Tz„  cosh  $£) 


sinh  C— ) 


sin  0. 


CA58) 


Substituting  these  expressions  into  the  horizontal  wave  force  equation 
yields: 


Fh  " CD 


pAn 


2T 


2H2  cosh2 


. .2  ,2nd.. 
sinh*  ( — 


cos  6 | cos  6 


-Si 


pV2n2H  cosh  ( Lj 
sinh  C~) 


sin  6 


(A59) 


or  Fh  - CD  FDh  cos  0 | cos  6|  - CM  F^  sin  6 
where 


(A60) 
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2U2  COSh2  (^-) 


sinh2  (-—) 


(A61) 


„ _ pV2n2H 

FMh  * r~ 
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cosh  (—■) 


sinh 


(A62) 


The  expressions  Fp^  and  Fj^  are  constant  for  a given  set  of  wave  and 
pipeline  conditions. 

The  horizontal  conponent  of  the  wave-induced  force  can  also  be 
written  as: 


Fh  “ ^2  cos  ® I cos® I + ®2  S2n  ® 


(A63) 


where 


A2  - CD  FDh  * !'2  CD  P A %ax  l>Wl 
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The  derivatives  of  this  expression  taken  with  respect  to  the  unknown  ' j 

parameters  kj  and  B2  and  set  equal  to  zero  give  the  following  equations:  J 


+ 2B~  sin  8.  cos  8.  |cos  8.  | 

4 X 1 1 


- 2 F . (8. ) cos  6.  I cos  8.  I 
oh  i 1 l 1 


= 0 


(A70) 
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oh 


( e.>r 


"5b: 


= 2 cos  | cos  6^|  sin  0.^ 


+ 2 B2  sin2  0^ 

- 2 F , (8.)  sin  8. 
oh  v i'  l 

= 0.  (A71) 


n 

Since  £ sin  0.  cos  0.  |cos  0. | * 0 for  an  even  number  of  equally 
i=l  1 1 1 

spaced  intervals  0.  summed  over  a complete  wave  cycle,  the  resulting 
summed  expressions  for  the  derivatives  set  equal  to  zero  are 

J (cos  0jL  | cos  8J)2  - ^ Foh  (6i)  cos  0i  |cos  0±  | = 0 (A72) 

n n 

B-  l sin2  0.  - I F , (0.)  sin  0.  = 0.  (A73) 

2 i=l  1 1=1  °h  1 


These  expressions  are  easily  solved  for  the  unknown  parameters  and 
B2,  yielding: 


A, 
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l F . (0. ) cos  0.  I cos  0.  I 
L . oh  v i ' i 1 i ' 
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l (cos  0.  (cos  0 |)2 
i=l  1 1 


(A74) 
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B - i = 1 
2 n 


7 F , (9.)  sin  6. 
L%  oh  ^ 1 i 


7 sin2  0. 
i=l 


The  coefficients  of  mass  and  drag  which  best  fit  the  horizontal  wave 
force  data  throughout  the  entire  wave  cycle  can  thus  be  obtained  directly 
from  the  parameters  and  since 


lD  = F, 


2 i=l 


it 

I F , (0. ) cos  6.  cos  6. 
. L.  oh  l i 


FDh  (cos  6i  lcos  eJ): 
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IF 
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PROGRAM  SMLOt A ( INPUT .OUTPUT » PUNCH) 

01  MENS  ION  X ( A | )tV(4U'2(4l),VS(41  > , V I ( 4 1 > .F  1(4 1 > ,FP1 4 1 > ,F< 4 1 1 . 
IF VI 4 I I ,RES(4l > ,0(131) ,P|4»  » ,0(4| I ,HI (41 ,MX(4) 

C SET  TEST  CONDITIONS 
ANGLE *0. 

0-2.000 

C READ  IN  DIGITIZED  DATA 

8 READ  I ,UP,DN,UF,OF,OIA,XC 
IF  I <JP  >10,10,9 

9 RE AO  2, CL ,T ,N, XV , XF ,C , ■ J ,FO , (HI  I 11 ,!■! ,41 
READ  3, (Ft (I) ,1-1,40) 

1 FORMAT(2F3. 3, 2F3.C, 2F3.3I 

2 FORMAT (F 3. 3,F3.2t  !2»2F 2-1 , 7F 3.0) 

3 FORMAT ( 24F  3*0 ) 

C DETERMINE  HAVE  HEIGHT 
DO  II  l>l|A 
IF ( H I ( I1-W0I2I ,23,23 
21  HX( I >--(VO-Hj ( I I )*( on/CIA 1 1 ./XV) 

GO  TO  11 

23  HX (l)a(HI(I) -VO ) A ( UP/C )*(l./XV) 

1 1 CONTI  N 3E 

H*>  (HX ( 1 ) ♦HX(2)-HX(  3 >-HX(  A)  ) /2. 

C CALCULATE  CONSTANTS  IN  FORCE  EOUATION 
PI* 3* I 413926636 
R-1.93A 

CALL  A A VEL ( T * D , XL 1 
ZV*CL ♦ ( » S*D I A I 
A-EXP( (2.*PI*ZV)/XL) 

COSHA*  A4 ( l./AI 
SINHA* A— ( la /A | 

8-EKPI (2* API 40) /XL) 

SlNHd^B-(  1./B1 
CSaCOOHA/SINHH 
SS-Sl NMA/SI NHB 

CFl-RAOl AAXCAh«PI AP I/(  2.ATAT) 

FLV-CF1*HACS*CS 
F0V«CF1 AMASSASS 
FMV»CFl*OIAAPI ASS 
U-(MACS*PI 1/T 

C DETERMINE  AVERAGE)  MAXIMUM)  AND  MINIMUM  FORCES 
FMAX  — .5 
FM IN**  S 
SF«0. 

DO  13  1-1,40 
IFIFKI  1-F0ir0,7l,71 
TO  FP(I»*(FIII >-FO)*(DP/C)4,002204/XF 
GO  TO  73 

71  FP(  I l*(FI( n-F0IA(UF/C>A,00220A/XF 
73  CONTINUE 

IF (FP( II *GT. FMAX >902)901 

902  FMAX-FP(I) 

901  IF(FP( I >.LT.FM!N1903)I2 

903  FMIN-FPI I 1 

12  SF-SF4FPII1 

13  CONTINUE 
SF-SF/AO* 


ri 

h 

} 


CA^CjfLATE  SUMS  OE  SOUARES  AMO  PRODUCTS 

SEX*Q? 

SEY-O. 

SEQ-O. 

SE2*0. 

OT a.  31  A l 5916936 
A*-DT 

00  IS  I • 1 1 AO 

A*A*OT 

XUt-COS  I A I 

««2>*A 

V«ll«CQS(St 

0(ti-stNiai 

C*S INI A I 


J(tl«C*»BS  1C) 
E«  I»*EP1I»-SE 


ipp*sFpmi)*Fm 

tPX*SPX»MU*X(ll 

ip»«»py*xM  »*y<  1 1 

SEQ-SFO*EI I )*Q<  1 1 

sf«spz»pin*im 

IS  continue 

AX-SEX/20. 

AYaSEV/20. 

AQmSEQ/20. 

A2-SEZ/I9. 

VX-AXASFX 

VY-AY*iEY 

VO-AQASEO 

VZ*AZ*SEZ 

CALC UCATEV*oe EPIC I ENTS  AND  PARAMETERS  EMI  ANO  K 
PHla28.64789*ATAN2I AO.AYI 
lEIPHI.LT. -AS. 

7999  PHl »PHI ♦ I EO. 

S999  CONTINUE 

» A • SORT  I A a*  AOA  AY*  AT  ) 

CLV-2. AVA/ELV 
ANG-(ANOLE*PI I/ISO. 

CLVA-CLV/COSI ANS I 
CLVU-CLVA/COSIAMGI 
CMV— — AX/FMV 
COV  — AZ/FOV 
XK«.5-ISE/ICLV*ELVII 
C PRINT  RESULTS  OE  ANALYSIS 
PRINT  200 
200  FORMAT 1 IHI I 

• i£&aS?8 

JOO^EORMAT  1 IS  \l  i 3.’ £e? 30^1  SIz  5E  i S^F  1 3.  I /// / » 

307  FORMAT  1 6X 1 1 2HT0T  AL  SUM  SO ,«X . 5MC0S2 A , » OX.SHSI N2A, IlX.AHCOST,  SX,  IOH 


xJJ.  !-'owt ',.;  JfWTJISfSS 


Kts.wite>r - 

301  F0*«TII*Ki»Ft6«*/l 


|H»ttl]X|IHir4l 


100  P0MAT141X,3MPLV,t2X«3HPMV,12X,3HP0V> 

mint  »c v«r»*v ••'ov 

303  «*»XT(J)X*JFIti»////»  • 

10*  PORNATI I0X*4H  CL V * l OX t St*  CLVA»10X»8H  CLVU, IIX*SH  CX»  ,I0X,8H  COV  t 

*RRI'5TH3©4,CLi»CLTA»CL¥U»CNT,eOV»XK»PHI 

308  FOTMT(SF»3.JiFl»«*iF«»tl////l 

300  PoiM4TiSOX|6HPAV6ILBI,TXf8rtPNAXILS»  » TX,6MPNINILBI» 

PUNT  304,SP.PMAX,PNtN 

3C4  KnCMT3ST?CL,OIaJaN6LB,T,m,XL,U,CLV.CLVA,CU|U,PHIjXK,CMV|COV 

PLGt’cRiIiNAL  OAT  A ANO  RESULTS  PCS  COMPARISON 

28  P0RMAT?*X,7M  FPU81|W»jh  P»(W)|M|*H  «C*(t#)F) 

OO  it  L— I • I 01 
G(kl*H 
31  CONTINUE 

PW1 l?*AX»it t 1*A0*0C 1 l*AV#Tt 1 l*A2*lt 1 >*SP 
RCS(  I I-PPIII-PVIII 
tS  CONTINUE 

e-AESIPPIIM 
OO  89  1-2,40 
C-ABSlPVI III 
A<ABS<PP<UI 
1P(C-B>87, 87,86 
S6 

»T  IP  I A-B 1 80, SO, 88 
88  8»A 


IP  188-BOO. 161,43,43 
81  1P1BB-400. 162,44,44 
8«  IP1B8-200. 163.48,48 
63'  IP <88- 1 00,  164,46,46 
64  88-80. 

60  TO  47 

43  88-800. 

60  TO  47 

44  88-400. 

60  TO  47 

48  88-200* 

60  TO  47 
46  88-100. 

*7  CONTINUE 

DO  32  1-1,40 
618II-IHI 
J.SI.»«8*PPtn 
K-St.*BB*PV<!> 

L-S1.t08*RES( 1 1 
6<J»— IH6 
6(K)«1H* 

pR,iNTl^lo|PP^^•2''^*^"*^<J,*• 

100  PORNAT 1 1HI,3P10.8, 101AI/I 
6IJI-IM 
6IKI-IH 
6ILIMH 
3*  CONTINUE 
60  TO  8 
10  CONTINUE 


Subroutine 
8-32. 4lS*3 1 " 

TPO-6. 2831 88*0 

JBHSU  INITIAL  ESTIMATE  POR  PAVELEN6TM 

IhaSop°RATER  INITIAL.ESTINSIE  POR  BAVELEMttTH 
3 *^_.T* SORT <0*32.  21 

* xl2b2tanm<  tpo/xlxi 

1PIAB8  1 XLX-XL I -.008 I 8, 4, 4 
8 RETURN 
END 
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<l_.« 


-HNINUk 


•01 

r>o 

iu 


•it 
• 12 
•OB 


Hl«t 

I*I*1 


m 


Mlll)iST,«.t«IO|l9l 
• 31  «l 


KIM 
*!■! 

IP<HI ( I I ,GT*0, >•!! .801 

iP(Ut-Tl(M3il,6t.iO,lat2,a32 

t*mi 

i»i*i 

1MHM  I j.GT.O,  1008, 113 


•13  T2(03l  _ _ 

60  TO  (814,8061, M3 
014  *3>M3M 

*I»,3*(xTl(2|-XTMn  > 

IX»K1 
I«IMX 
60  TO  901 

*06  TM1>( T2(1)+TI( II 1/2* 

TM2>(T2(2I+TI(2) 1/2, 

(TNKT9IM 

XT*2*TN24»6 

IMKMW1 

ITM2-XTM2 

TT1MTM1 

tT2>ITN2 

T-(TT2-TrtJ4,02 

X/oflO,*T 

UKlJld 

J>XJJ 

CALCULATE  C9N0T4HTS  IN  P0*CE_  EQUATION 
PI >3,1419926936 
*>1,938 

call  wavel(T,d,xli 
ZV,CLM.S401A> 

A>EX»( (2,4Pt*ZV)/XL» 

COSHA>A*(l./AI 

SIKHA>A-<1./4I 

••>BXP((2,*PI*9|/XL) 

S1NHS>9-(|,/S) 

CS>COSHA/SINHB 
S5>S INHA/S I NH8 

CPt>ROt>(**XC6H4PI*P|/(2,*T*T> 

PLV-CPt*H*C9«CS 
POV>CP|4H45S*SS 
PM V»CP 1 40 1 A 4P I 458 
W-(H*C5*P! >/T 

DETERMINE  AVERAGE,  MAXIMUM,  AND  MINIMUM  PCRCE3 
PMAX  — ,3 
PMIN>,5 
SP>0, 

00  1 3 I>t,J 

PP( I »>P| (IaITMI- 11*0,0220* 
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m 


ttT*****  **•*  •** 1 


t > 

CONTI NUt 

| ftt?  *'"*  00  MUMI*  otto  «fl#«CT* 

Si:?: 

•oo*o. 

ot is? i ♦ i ‘ * 
a«-ot 

00  IS  !-»»■> 

IMUtDI 

Xttt*C0S  1*1 

VI I l*C0SIS> 

•|lf«StN(S» 

... 


IC» 

JWWJUlTft.... 
8SSS?}»«!!I 
8B8S8S!!fini 
**  a;::«*o 

K5::s::i , 

H£»J3KJ 

AV-SOV/SW 

MfVI/UI 

sasua 

TSSS 

***'  yS**oS*T^  to* AO**V**V  • 

CC*A«CI.V/CO»IAN«l 

CNV—AX/OOV 

Sirsif  so/I O.VSW.V » > 


Of  MMiTtlt 


m? 

Ml 

III 


4R^ipWMr~ 

winrTMi  ' 1 W 

MJ  rOMMTIM**  Ml ••«////» 

MINT  194 

4M  MMITIIMilH  aV|IU«lN  CU*A»it *t«M  <XVU«ttX*tH  CM*  »IM||N  COW  « 



jit  ooSSat  (Son. omo a toils  » « t«  ,oho ha  m l* i •rx»0M9*t««tL9> » 

MINT  104.il*  ,W44*.94JN 
J04  MNMITI  Jlhl Iltl (ii|l«l//// 1 

rUHCH  MT»CL»OtA«AMSLKt  T »H*  «.*  Uf  CUY.CLVA.CLVM.PMI .M  .CHW.CO  V 


C KCT  MUtlW.  OAT*  AMO  MIIATI  90*  COMMIHH 

~«*i  wiTtirlitJliiiitiS' 

Ml  Jt«S 

*0  TO  904 
AOS  JtO 

SO  TO  99* 

Ml  JIM 
Ml  CaMMMMI 
-MMMFNIHI 
04AN4XKA.CI 
M>».  /* 

CClllt/NNIXll) 

MINT  M 

14  rOIMTI«X,TH  Ma»,MiTN  FWall»l«tIN  Rtllkll/I 
00  91  L*l,lll 
OILI-tM 
JU  eONTIMUC 

00  30  l*Wi|Jl 


OtO»>*lMl 


•Mil  1 4IT01-1 I 
Ml  I )>AX*Ktt  lTA0*0il  14AV9TC  IHimil  IHF 
M«UI#F(II*MI  I* 

M*ll.*M«F»(ll 


4I4JI" IH* 
IIMlwlH* 

ML 


.o.  Si5iT!Wi:y.i::r.iw5ru^ 

mm 

•ILL  |*1H 
0(011  «IH 

*wyrf 

to  COM ft MOO 

0^r*»AtS^  IhUIm,  «MlSATJt  900  NAySUWATH 

* farawffttt........ 

•»  BT" 
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I 


PRUORAM  H0MI2L t I NPUV, OUTPUT .PUNCH; 

Ot«*CNSIOM2<  All  .PHI  >,YS<*1  > *Yl  141  >,P1UI>,FP|  41  l.FIAt t .FHI41I. 
IMSItl  I iCI  tOtt 
set  test  ccnoiticns 
ou«c«jii 
XCOaSir 

ANGLE *0* 
i>«2»C00 

A1AO  IN  0IGITI2CD  OATA 
S ME  AO  1. UP. ON. CFO, CPU 

ifiupI 10,10.4 

* MAO  2,CLfT,N*X«tXF,C,M0,FO 
MAO  3.IYI  < I I .XI  « t I ,1*1  ,401 


1 FORMAT < 1F3. 3.2F 3.0) 

« FORMAT <*3*  3 .P3»t , IS . 2F1. I , 3F3,0 1 
3 roMAT(i*F3.0) 

C DETERMINE  RAVE  HEIGHT 
00  II  1>1, 31.10 
IFtYtl I )-VO)21, 23.13 
it  ys(ii«*(M>«t(ii)iioN/cim>/xii 
CO  TO  It 

13  TSIU4rl(l)-lOI*IUP/CI*(U/MI 
11  CONTINUE 

h.<ys< ii»vs<ti>-vs(t i>-vs<  st>  >/2» 

C CALCULATE  CONSTANTS  IN  FORCE  EQUATION 
PI  *3* 1419926834 
MI.13S 

CALL  PAVELI T.O.XL) 

1V*CLM,0*OIA> 

AaEXPi <£«*PI*2Yl/XL> 

CQSMA-A«< I* /A I 
R»EXP< <1,*PI*0|/XL» 

SINHO-B-<l./*> 

CS-CaSHA/StNNB  

CP|>A*0IA4xC*H»Pt*P|/<*.«T*T) 

fon«cfi*h*cs*cs 

PNH-CP 1 *0 I A*P t *CS 
U«(H*CS*PI l/T 

C OETERM INC  AVERAGE*  MAXIMUM*  ANO  MINIMUM  FORCES 

FMAX  — .3 
FNIM-.S 
SF-Q. 

00  13  l-l.AO 
IFIFK  I 1— FO  ) TO  1 . TOC* T02 
rot  FPin-iFo-Fin  »Ip<c#o/c»*.ooiioaxxf 
GO  TO  T03 

T02  FPU  >■-<*! (tl»FOI«<CFU/CI**Q02»04/XP 
T03  CONTINUE 

IFlFPt I >.GT.FMAX)902,40I 
ROE  FRAX«FPm 
*01  IFfFPf I I. LT.PMIn>903,12 
*03  FMIN-FOtn 
IS  SF-SF*FP«I» 

13  CONTINUE 
SFaSF/40, 

C CALCULATE  SUMS  OF  SQUARES  AND  PROOUCTS 
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SFFao. 

SFP-O. 

SFZ-O. 

0T«. 31*13926536 

»■- or 

00  15  I»1,*0 

»«**0T 

PI  I )>S*N(AI 

CaCOSIAl 

ZII)-C*AB3  (C) 

FI  I I aFPI t » 

SFFaSFFaF { I >*F|I ) 

SFP>SFO»F| I >*P| I ) 

SF2aSFZ*F| I >*Zi I) 

IS  CONTINUE 
APaSFP/20* 

AZaSPZ/lSa 

MP»AP*SFP 

VZaAZ*SFZ 

CALCULATE  COEFFICIENTS 
CMM-»-AP/FMH 
COM» AZ/FOH 

PRINT  RESULTS  OF  ANALYSIS 
PRINT  200 
200  FORMAT! IH1) 

PRINT  « . ... 

6 FORMAT! 1 OX, 6HTI SEC) , TX , 6HHT I F T ) ,SX ,9HBA VEL IFT I *6X • 7H0EPIFT I , AX, 9HU 
1 MAXI FPS) ,4X,9HCLEARIFT) , 6X, 7HO I A< FT ) , 3X , I2HCVL  LGTMI FT  I , AX, BHANG! O 
2IC)  I 

PRINT  300,  7. H,  XL  ,D,U»CL,QIA»XC*ANGLE 
300  FORMAT I2X.F13#2,FI3*3,F1 3,2 , SF 13, 3,F 1 3a I////I 
PRINT  330 

330  FORMAT  I 26X, 12HTOTAL  SUM  S0,7X, 3HSIN , I0X, OHCOS/COS/ 1 
PRINT  331  LSFF.YP.YZ 

331  FORMAT I 20X , 3F ISao/ ) 

PRINT  3 32 

332  FORMAT I 46X, 2MAP , 13X ,2HAZ I 
PRINT  333.AP.AZ 

333  FORMAT  4 3ftX , 2FI 5a  8/ • 

PRINT  33* 

33*  FORMAT |46X , 3HFMH , 12X , 3HFDH) 

PRINT  JSSL'XNiFOH  . . . 

335  FORMAT  1 3SX, 2FlSa  6//// 1 
PRINT  336 

336  FORMAT I46X, 3HCMH, 12X ,3HCDHI 
PRINT  337,CMH,CQH 

337  FORMAT  1 35X , 2FI Sa  J/ /// I 
PRINT  309 

309  FORMAT I3SX , SHFAVGI LSI , 7X , BHFMAXI L8 I , 7X, 8MFMINILS I > 

PRINT  304,SF,FNAX.PM1N 
304  FORMAT! 30X,F15.6,2F|5.S////) 

PUNCH  367, CL, 01  A , ANGLE, T , H , XL ,U , CMH.COM, SF 
387  FORMAT (F4. 3 , F5, 3,F4,0,F9,i,FS.3,F6*2,F6, A, 2F3,  2,FI0*6) 

PLOT  ORIGINAL  OAT  A AND  RESULTS  FOR  COMPARISON 
PRINT  26 

26  FCRMAT|4X,7H  FPILB),3X,7M  FHILB),2X,SH  PKSI'.S)/) 


1 


AaABStPMIN) 

s:au;?Ki, 

BBa JS./8 
DO  31  L«1,10I 

31  CONTINUE 

DO  32  !■!.*<> 

JHt  t |aCOH*FOH*2( I )-CMH*FNH*P< I ) 

Rcs<n«FPm-PH«i> 

6(St»a|HI 
J*5 1 • *8B*PP ( 1 J 
K>$Ii*I8*FH( t I 

6(X|a|H4 

CILtalH* 

..o  aawtwa'.ijRSisiBr ; ■»•*- 

G( J|a|H 
G(  K ) a 1H 
6<L1«1H 

32  CONTINUE 
40  TO  B 

10  COMTINUE 

EMO 

SUBROUTINE  NAVEL ( 0,XL) 

B«  32* I*T*T/». 283143 
TPO-6. 283185*0 
If IB-TPOl  2,2.31 

OeEP  MATER  INITIAL  ESTIMATE  VOR  WAVELENGTH 

2 XL  “8 

60  TO  4 

8HAU.Q*  EATER  INITIAL.  ESTIMATE  fOR  WAVELENGTH 

3 XL • TAJORT I 0*32, 2 ) 

4 XLX-XL 
XL-B4TANNITPO/XLX) 

If  IABS  IXLX-X'-  >-.005)8,4,* 

5 RETURN 
ENO 


i 

) 
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APPENDIX  E 

TABULATED  VERTICAL  FORCE  DATA 
FROM  TWO-DIMENSIONAL  EXPERIMENTS 


t- 

ty'v ' 


i: 

4. 


V ‘ \ 

: 
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cur 

Dl  4 

4MG 

T 

H 

L 

UMAX 

CL  V 

CLV* 

cl  vu 

PHI 

K 

CM  V 

cov 

a AO  1 

<UJ 

0. 

1*13 

• ooe 

7.26 

• 0427 

0*  14 

0,  14 

0,14 

26.20 

• 4850 

2*60 

-818*75 

• ooa 

• 333 

0« 

1.23 

• 117 

7*24 

• 110? 

7,  60 

7,60 

7.40 

29.43 

.0974 

2.99 

* 278*61 

• COI 

• 333 

0. 

1 • 21 

• 147 

7.14 

• 1554 

v,6t 

7.02 

7.02 

24.87 

•0456 

<•74 

* 169*  30 

•cot 

• 333 

0. 

1.24 

• 207 

7.37 

,1990 

4.  26 

4*24 

0.  26 

10.67 

• 5428 

1.4* 

-168,70 

• 001 

• 333 

Oi 

1*24 

• <57 

7,37 

• 2471 

4,03 

4,5  3 

4.53 

19*74 

• 8064 

2*91 

-117,97 

• 001 

• 333 

0* 

1 • 24 

.27? 

7.37 

• 2647 

4.62 

6,02 

6.6  2 

10.27 

• 6174 

2*49 

-114,63 

• 001 

• 333 

0. 

1 • 2 A 

• 295 

7.57 

• 2932 

6,42 

4*42 

6.42 

21.01 

• 4047 

3*15 

-03,52 

• cot 

• 333 

0. 

1.45 

• 044 

9.  39 

• 0743 

3*44 

3.C  4 

1,64 

22,47 

• 0122 

2,99 

-202*49 

• 001 

• 333 

0. 

1.45 

• 109 

9*39 

• 1341 

4*01 

0*01 

4.01 

21.17 

• 4970 

2,44 

-127,40 

• 001 

• 33  J 

0* 

l *4  A 

• 144 

9.44 

• 2029 

5,  74 

8*76 

5.76 

ta.09 

• 4444 

<•70 

-64,73 

• 001 

• 333 

0« 

U49 

• 101 

9.39 

• 2500 

5*  14 

8*  14 

5.14 

14.49 

• 4017 

1.40 

-79*81 

• 001 

• 33  3 

0* 

1.4  f 

• <04 

9.50 

• 2077 

4,94 

4,94 

4.94 

14*63 

• 3972 

3*29 

-91,25 

• 001 

• 333 

0* 

1.4  7 

• <34 

9.50 

• 2901 

5.  •< 

6.92 

5.92 

13.50 

• 3894 

1*44 

-64* 12 

• 001 

« 33  J 

0* 

1*46 

• 244 

9.67 

• 33*0 

5.35 

5*35 

5.35 

13.  13 

• 3422 

3.20 

-39,00 

• 001 

• 333 

0. 

l • 4 A 

• 249 

9.40 

• >324 

5*44 

5*04 

5.04 

10,98 

• 3817 

2*  14 

-96*40 

• 001 

,333 

0« 

1*49 

• <03 

9*77 

• 3594 

5*75 

5,75 

5.75 

14.97 

• 3030 

2*90 

-32*33 

• 001 

• 333 

0. 

1*02 

• 051 

1 2*  AO 

• 0774 

4,40 

4*08 

4,05 

23*49 

•8744 

3*45 

-730.11 

• 001 

• 333 

0* 

1.02 

• ooo 

12.00 

• 1330 

8*40 

8*00 

5,00 

19*37 

*4447 

4*23 

-444,46 

• 001 

• 333 

0. 

1*02 

• 007 

Iff. 40 

• 1319 

4,17 

4.17 

4.  17 

20,09 

• 9349 

2*49 

-410.  12 

• 001 

• 333 

0* 

i»at 

• 127 

12,7| 

• 1910 

4*04 

4*04 

4,04 

13«70 

• 4214 

2*74 

-170.I9 

•not 

•333 

0. 

1,52 

• 145 

12.10 

• <405 

4*14 

4,  18 

4.11 

16*27 

16478 

2*87 

-270*04 

• 001 

• 333 

0. 

1.05 

• 190 

13,07 

• <094 

5,  41 

8«4l 

S.4t 

14.61 

• 3324 

3.4  2 

-212,08 

• 001 

• 333 

0* 

1.05 

• 217 

13,07 

• 3317 

8*94 

8,98 

6.95 

12*02 

• 3437 

<•21 

-204,74 

• 001 

• 333 

0* 

1.06 

• 241 

13*07 

• 3403 

8,  88 

5,55 
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APPENDIX  g 

TABULATED  HORIZONTAL  FORCE  DATA 
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